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Abstract 



The evolution of the usual Erdos-Renyi random graph model on n vertices can be described as 
follows: At time start with the empty graph, with n vertices and no edges. Now at each time k, 
choose 2 vertices uniformly at random and attach an edge between these two vertices. Let G„(fc) be 
the graph obtained at step k. Following |22j and [25], the work of Aldous in [4] shows that for any 
psj I fixed t S R, when k{n) = n/2 + the sizes of the components in G„(fc(n)) scale like n^/'^ and, 

^ ■ as n — > oo, the scaled component size vector converges to the standard multiplicative coalescent at 

(N ■ time t. Furthermore, when k{n) — sn/2, for any s < 1, the size of the largest component is O(logn) 

^^i^ ■ while for s > 1, the size of the largest component scales like f{s)n, where /(s) > 0. Thus = 1 is the 

critical parameter for this model at which time a giant component emerges. 

The last decade has seen variants of this process introduced, loosely under the name Achlioptas 
processes, to understand the effect of simple changes in the edge formation scheme on the emergence 
of the giant component. Stimulated by a question of Achlioptas, one of the simplest and most popular 
of such models is the Bohman-Frieze (BF) model wherein at each stage fc, 2 edges ei(fc) = (wi,W2) and 
62 (fc) = (^3,^4) are chosen uniformly at random. If at this time vi, V2 are both isolated then this edge 
is added, otherwise 62 is added. Then [T2^ (and further analysis in IT) show that once again there is a 
. critical parameter, which is larger than 1, above and below which the asymptotic behavior is as in the 

Erdos-Renyi setting, in particular this simple modification in the attachment scheme delays the time 
of emergence of the giant component. While an intense study for this and related models seems to 
suggest that at criticality, this model should be in the same universality class as the original Erdos-Renyi 
process, a precise mathematical treatment of the dynamics in the critical window has to date escaped 
analysis. In this work we study the component structure of the BF model in the critical window and 
show that at criticality the sizes of components properly rescaled and re-centered converge to the same 
limits as for the original Erdos-Renyi process. Our proofs rely on a series of careful approximations of 
the original random graph model leading eventually to the treatment of certain near critical multitype 
branching processes with general type spaces through functional and stochastic analytic tools. 

Key words: critical random graphs, multiplicative coalescent, entrance boundary, giant component, 
branching processes, inhomogeneous random graphs. 
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1 Introduction 



Random graph models of various systems in the real world have witnessed a tremendous growth in the 
last decade. The availability of a large amount of empirical data on real world networks and systems such 
as road and rail networks, bio-chemical networks, social networks and data transmission networks such as 
the internet has stimulated an inter-disciplinary effort to formulate models to understand such systems 
ranging from (largely) static systems such as road and rail networks, to more dynamic systems such as 
social networks, information networks and models of coagulation and aggregation in physics and colloidal 
chemistry. 

The classical Erd6s-Renyi(ER) random graph can be thought of as a network evolving in time via the 
following prescription: Start with the empty graph at time zero. Then at each discrete time point, choose 
an edge uniformly at random and place it in the network (irrespective of whether it was present before or 
not). This simple model exhibits a tremendous amount of complexity, see fl5ll27j . As one adds more and 
more edges, eventually the system transitions from a "sub-critical" regime wherein the largest component 
is of order logn, to the super-critical regime wherein there exists a unique giant component of order n 
and the next largest component is O(logn). Understanding the emergence of this giant component and 
the properties of near critical Erdos-Renyi random graph have stimulated an enormous amount of work, 
see e.g [35l [HI [19] and the references therein. 

This model has been modified in various ways the last few years, to understand the effect of choice in 
the emergence of the giant component. For example, in p] and pOj simulation based studies are carried 
out for an attachment scheme where at each stage one chooses two edges uniformly at random, and then 
one uses the edge which delays the emergence of a giant component (for example by choosing the edge 
that minimizes the product of the two adjoining components). A number of fascinating conjectures about 
the "explosive emergence" of the giant component are stated in these articles. 

Such modifications of the ER random graph model have generated a tremendous amount of interest in 
many different communities, ranging from combinatorics to statistical physics. Although a number of such 
models have been analyzed in the sub-critical and super-critical regime, see for example the marvelous and 
comprehensive [H], understanding how the giant component emerges, even for very simple modifications 
to the ER setting, for example for the Bohman-Frieze model (which we describe below) has escaped a 
mathematical understanding. The aim of this paper is to prove the first rigorous results about such 
models at criticality and obtain precise asymptotics for the merging dynamics in the critical window that 
lead to the emergence of the giant component. 

The Bohman-Frieze (BE) model is as follows: Let 0„ be the empty graph on the vertex set [n] 
{1,2, ...,n} with no edges. The random graph process {G^^(/c)}fcgN = {G„,(A;)}fegN evolves as follows: 

• At time fc = let 0^(0) = 0„, 

• The process evolves in discrete steps. For /c > let the state of the graph at time k be Gn{k). Then 
at time A; -|- 1, the graph G„(fc -|- 1) is constructed as follows: Choose two edges ei = (fi,f2) and 
^2 = (^^3)^4) uniformly at random amongst all the possible (2) edges. If vi,V2 are isolated vertices 
then let G„(A; + 1) be the graph G„(/c) with edge ei added else let G„,(/c -|- 1) be the graph G„(A;) 
with the edge 62 added. 

Note that the standard Erdos-Renyi (ER) process can be thought of as a variant of the above process 
wherein we add edge ei at each stage (irrespective of whether the end points are isolated vertices or not). 
We shall use G^^(/c) to denote the usual Erdos-Renyi process. 

For the standard Erdos-Renyi process classical results (see e.g. [22], [21], [15]) tell us that by scaling 
time by n and looking at the process G^^([?ii/2j), a phase transition occurs at tc{er) = 1, namely 

(a) Subcritical regime: For fixed t < 1, the largest component in G^^([nt/2j) is O(logn). 

(b) Supercritical regime: For fixed t > 1, the largest component has size ~ f{t)n for some positive 
function / (namely there is a giant component) while the second largest component is of size O(logn). 

(c) Critical regime: For i = 1, the first and the second largest components both are G(n^/^). 
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For the BF model, the original work in p!^ shows that there exists to > 1 such that the size of the 
largest component in G^^( [nto/2j ) is Op{n). Thus this model "avoids" the giant component for more 
time than the original Erdos-Renyi model. The comprehensive paper [31] showed for the BF model and 
for a number of other related processes, that there exists a critical value tc, such that conclusions (a) 
and (b) hold for these models as well with 1 replaced by t^. For the BF model, the supercritical regime, 
namely the regime t > tc = tc{hf) was studied further in [26j. Numerical calculations from this paper 
show that tc{bf) ~ 1.17. 

Understanding what happens in the critical regime for these models is of tremendous interest for 
many different fields, as it turns out that much of the "action" happens at criticality. For the Erdos- 
Renyi processes, Aldous's result in [4j shows that if we denote the component size vector (components 
arranged in decreasing order), in the graph G^^([nt/2j), by 

Cr(t) = (C«(t):i>l), 

then, for any fixed A G M, as n — t- oo, 

Cr (A) n-^/^CT {^ + ^)^ (1.1) 

where the process {X{X) : -co < A < oo) is a Markov process called the eternal standard multiplicative 
coalescent on the space 

if = {{xi,X2,...) ■■ xi>X2> ■■■ >0,^Xi < oo}. (1.2) 

i 

The space /| is endowed with the same topology as inherited from and the above — ^ denotes weak 
convergence in this space. In fact the paper [3] shows weak convergence of C^^ to X in P((— oo, oo), /|) 
(the space of RCLL functions from (—00,00) to endowed with the usual Skorohod topology). We defer 
a full definition of the above Markov process to Section 12.11 

Despite the tremendous amount of interest in the critical regime of such models and the significant 
progress in analyzing these models above and below criticality (see e.g. [31] for a general treatment of 
such models and [26] for an analysis of the BF model), analysis of the component sizes at criticality and 
an understanding of the dynamic behavior of these components and how they merge through the critical 
scaling window for the giant component to emerge, for settings beyond the classical Erdos-Renyi model 
has remained an extremely challenging problem. The most commonly used technique in getting refined 
results about the sizes of the largest component at criticality, which is to show that the breadth-first 
exploration processes of components converge to a certain 'inhomogeneous Brownian motion' does not 
seem to extend to such settings and one needs to develop a rather different machinery. In this work 
we take a significant step towards the understanding of the behavior of such models by treating the 
asymptotics in the critical regime for the Bohman- Frieze process. Although this model makes a simple 
modification to the attachment scheme for the Erdos-Renyi setting, the analysis requires many new 
constructions and mathematical ideas. Our main result. Theorem 12. H shows that at criticality the sizes 
of components properly rescaled and re-centered converge, as for the original Erdos-Renyi process, to 
the eternal standard multiplicative coalescent. The proof proceeds through a sequence of approximations 
that reduce the study to that for a certain inhomogeneous random graph model which allows us to bring 
to bear techniques from the theory of multitype branching processes on general type spaces and integral 
operators on Hilbert spaces. These culminate in an estimate, obtained in Proposition 15. H on the size of 
the largest component through the "subcritical" window. Using this estimate and other careful estimates 
on quadratic variations of certain martingales associated with the Bohman-Frieze process we then obtain 
(see Proposition 14. ip precise asymptotics on the sum of squares and cubes of component sizes just below 
the critical window. The latter result is a key ingredient to the proof as it allows us to show that just 
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before the critical window, the component sizes satisfy the regularity conditions required for convergence 
to the standard multiplicative coalescent, and thus allow us to couple the BF process with a standard 
multiplicative coalescent through the critical window to prove the main result. Although not pursued 
here, the techniques developed in this work for establishing the asymptotic behavior in Proposition 14.11 
are potentially of use for the analysis of the critical regime for a very general class of Achlioptas processes 
(called bounded size rules). 

Organization of the paper: Due to the length of the paper, let us now briefly guide the reader 
through rest of this work. We begin in Section [2] with the construction of the continuous time version of 
the BF-process and then state the main result (Theorem 12. ip in Section [2.11 Next, in Section [3] we give 
a wide ranging discussion of related work and the relevance of our results. We shall then start with the 
proof of the main result by giving an intuitive sketch in Section [4] wherein we also provide details on the 
organization of the various steps in the proof carried out in Sections 5-7. Finally Section [8] combines all 
these ingredients to complete the proof. 



2 The Bohman-Frieze process 

Denote the vertex set by [n] = {1, 2, . . . , n} and the edge set by <f„ = {{vi, ^2} ■ vi ^ V2 £ [n]}- To sim- 
plify notation we shall suppress n in the notation unless required. Denote by BF(t) = BF„(t), t £ [0, 00), 
the continuous time Bohman-Frieze random graph process, constructed as follows: 
Let S"^ = £ X £ he the set of all ordered pairs of edges. For every ordered pair of edges e = (ei, 62) G 
let Ve be a Poisson process on [0, 00) with rate 2/n^, and let these processes be independent as e ranges 
over £'^. We order the points generated by all the (2) x (2) Poisson processes by their natural order as 
< ti < t2 < •••• Then we can define the BF-process iteratively as follows: 

(a) When t £ [0,ti), BF(t) = 0„, the empty graph with n vertices; 

(b) Consider t E [tk,tk+i), k £ N+, where is a point in "Pg and e = (61,62) = ({^^i, ^2}; {"^3; ^4})- If vi, 
V2 are both singletons (i.e. not connected to any other vertex) in BF(tfc— ), then BF(t) = BF(tfc— )U{6i}, 
else let BF(t) = BF(tfe-) U {62}. 

Note that multiple edges are allowed between two given vertices, however this has no significance in our 
analysis which is primarily concerned with component sizes. 

Consider the same construction but with the modification that we always add 61 to the graph and 
disregard the second edge 62. Note that the total rate of adding new edges is: 

n\ / n\ 2 n 

2) V2j;^^ 2- 

Then this random graph process is just a continuous time version of the standard Erdos-Renyi process 
and the convergence in (jl.ip continues to hold, where C^^(t) now represents the component size vector 
at time t for this continuous time ER process and once more, tc = 1 is the critical parameter for the 
model. 

As proved in [41], the Bohman-Frieze model also displays a phase transition and the critical time 
tc ~ 1.1763. We now summarize some results from the latter paper that characterize this critical parameter 
in terms of the behavior of certain differential equations. 

The following notations and definitions mostly follow [26j. Let Cn\t) denote the size of the largest 
component in BF„(t), and Cn{t) = {Cn\t) : i > 1) the component size vector. For convenience, we define 
C4'' (t) = whenever t < 0. 

For fixed time t > 0, let Xn{t) denote the number of singletons at this time and x{t) = Xn(t)/n denote 
the density of singletons. For simplicity, we have suppressed the dependence on n in the notation. For 
A; = 2, 3, let 

5.(i) = E(^n'(*))' (2.1) 

i>l 
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and let Sk{t) = Sk{t)/n. Then from [3T], there exist deterministic functions x[t) , S2{t) ^ s-i{t) such that for 
each fixed i > 0: 

x{t) — )■ x{t), Sk{t) — > Sk{t) for k = 2,3, 

as n — )• oo. The Hmiting function x{t) is continuous and differentiable for ah t G M+. For k > 2, there 
exists 1 < tc < oo such that Sk{t) is finite, continuous and differentiable for < t < tc, and Sk{t) = oo for 
t > tc- Furthermore, x, 82,83 solve the following differential equations. 



x'{t) 


= -x'^{t)-{l-x'^{t))x{t) 


for t G [0,00,) 


x{0) = 1 


(2.2) 




= X\t) + {1 - X\t))8l{t) 


for t £ [0,tc), 


S2(0) = 1 


(2.3) 




= 3x^{t) + 3{1 - x'^{t))s2{t)s3{t) 


for t £ [0,tc), 


53(0) = 1. 


(2.4) 



This constant tc = tc{bf) is the critical time such that whp, for t < the size of the largest component 
in BF(t) is O(logn), while for t > tc there exists a giant component of size 0(n) in BF(t). Furthermore 
from [26] (Theorem 3.2) there exist constants 



a = (1- x'^{tc)y^ ^1M3 (2.5) 



/3 = lim-^^Ri.764 (2.6) 



such that as 1 1 



S2{t) ~ (2.7) 



tc-t 

\3 



a3 



2.1 Main results 



duced in Section [TJ For x G let ord(3;) £ if he a reordering of the vector x in the decreasing order. For 
X G /?, 1 < i < j < 00, let 



Our goal in this work is to establish a limit theorem of the form in (jl.ip for the BF process. 

We begin with a precise definition of the standard eternal multiplicative coalescent process X intro- 

x^^ = ord(x + {x,ej){ei — e^)), 

where {ej} is the canonical basis in P. Namely, x*-' is the vector obtained by merging the i-th and j-th 
'components' in x and reordering the vector. Aldous(j4]) showed that there is a Feller Markov process 



with sample paths in D([0, 00) : if) with infinitesimal generator 



Aucfix) = Y,XiXj{f{x'n - fix)), xGlf, f -.if^R. (2.9) 

i<j 

This Markov process describes a coalescence dynamics where at any time instant the i-th and j-th clusters 
merge at rate equal to the product of the sizes of the two clusters. One special choice of initial distribution 
for this Markov process is particularly relevant for the study of asymptotics of random graph models. We 
now describe this distribution. Let {W{t)}t>o be a standard Brownian motion, and for a fixed A G M, 
define 

Wx{t) = W{t) + Xt--, t > 0. 
Let Wx denote the refiected version of Wx, i.e., 

Wx{t) = Wx{t) - min VF'a(s), t > 0. (2.10) 

0<s<t 
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Define an excursion of Wx as an interval {l,u) C [0, +00) such that Wx{l) = Wx{u) = and Wx{t) > 
for ah t G {l,u). Define u — I as the size of the excursion. Order the sizes of excursions of Wx as 

ei(A)>e2(A)>e3(A)>--- 

and write H(A) = (Cj(A) : i > 1). Then H(0) defines a Z| valued random variable. Denote by {X(A)}a>o 
the /| valued Markov process with initial distribution as the probability law of H(0) and infinitesimal 
generator as in (j2.9p . Then [3] shows that for each A € (0,oo), X{X) has the same law as H(A). In fact, 
[5j shows that the process X can be extended for A G (—00,00), namely there is a process with sample 
paths in P((— 00,00) : Z|), denoted once more as X, such that for every A G (—00,00), X{\) =d S(A) 
and {X{\ + t)}t>Q is a Markov process with generator ^mc- This description uniquely characterizes a 
probability measure on P((— 00,00) : /|) representing the probability law of X. The process X is called 
the eternal standard multiplicative coalescent. 

We are now in a position to state our main result. For the rest of this work tc = tc{bf) will denote 
the critical point for the continuous time Bohman-Frieze process as defined above. Recall the constants 
a and (3 defined in (US} and I^M). 



Theorem 2.1 For XgR, let 



(A) = [t. + /5^/3a-^ j : ^ > 1 j (2.11) 

be the rescaled component sizes of the Bohman-Frieze process at time tc + af3'^^^^^^Y/3- Then 

Cf,^ X 

as n —7- 00, where X is the eternal standard multiplicative coalescent, and — ^ denotes weak convergence 
in the space P((— 00,00) : /|). In particular for each fixed A G M, the rescaled component sizes in the 

critical window satisfy C'^^(A) — ^ '='(A) where H(A) are the sizes of the excursions for the reflected 
inhomogeneous Brownian motion in ()2.10p . 



3 Discussion 

Let us now give a wide ranging discussion of the implications of the above result and proof techniques. 

(a) Dynamic network models: The formulation and study of dynamic networks models, such as the 
Bohman-Frieze process wherein one has both randomness and choice, are relatively new. Motivated 
by a question by Achlioptas, that whether simple modifications of the Erdos-Renyi random graph 
model could allow one to avoid the giant component for a larger amount of time, such models now 
fall loosely under the term Achlioptas processes (see e.g. [I2l [I3l HH [H]). The models considered 
in the above papers use bounded size rules, wherein there is a fixed K such that when making 
a choice between the two randomly selected edges, all component sizes greater than some K are 
treated identically (e.g. K = 1 for the model considered in this study). Much recent interest has 
been generated by recent models formulated in [Ij, wherein interesting examples of unbounded size 
rules have been analyzed. Simulations of such models seem to suggest a new phenomenon called 
"explosive percolation", wherein the phase transition from a subcritical to a supercritical regime 
with a giant component appears much more "abruptly", in that the largest component seems to 
increase from a size smaller than ^/n to a size larger than n/2 in a very small scaling window. 
Examples of such rules include the product rule wherein one chooses two edges at random and 
connects the edge that minimizes the product of the component sizes on the 2 end points of the 
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edge. See [20[|18| for additional models of this type. Recently in f30] it was shown that for a number 
of such models, the phase transition is actually "continuous" , so for such models one expects similar 
behavior as what one sees in the Erdos-Renyi random graph model. However, understanding what 
happens at criticality and the behavior of the scaling window is a challenging math program. The 
techniques in this paper have the potential to be extended to the analysis of a number of such 
models, which we attempt to do in work in progress. 

(b) Recent results on the Bohman-Frieze model: After this work was submitted for publication, 
we came across the interesting preprint [29j that was announced around the same time. The latter 
paper studies the Bohman-Frieze process in the cases when t = tc — £ and t = tc + £, for fixed 
e > 0. In [29] the largest and second largest components are studied, bounds on the sizes of these 
components, as well as the number of surplus edges amongst all components in these regimes are 
derived. The techniques used in |29j for understanding the size of the largest component in the 
subcritical regime are very different from those used in the current paper. In particular, our work 
requires an understanding of the fine scale asymptotics of the entire vector of component sizes at 
criticality (properly rescaled). For identifying the scaling window and how the giant component 
emerges via the merger of small components through the scaling window, we need more refined 
results at the critical value, in particular we need to study the behavior as e = e(n) — t- reasonably 
quickly - see Proposition 14.11 where precise asymptotic results for e = 7 G (1/6,1/5) are 
obtained. Conjecture 1 in [29] states that one expects an upper bound of log n on the largest 
component at the time instant (tc — e). This should be compared with the (logn)^/e^ upper bound 
for times tc — e established in Proposition 15.11 of the current work, not just for a fixed e but for 
e = e{n) — t- 0. In fact, the proposition establishes an estimate that is uniform over the time interval 
(0,tc — ''^~^')^ 7 G (0,1/5). This proposition is at the heart of our analysis and its proof requires 
significant work and new ideas and techniques for general random graph models with immigrating 
vertices and near critical multi-type branching processes with general state spaces, see Sections [5] 
and[6l 

(c) Multiplicative coalescent: As shown in the current work and several other papers, see e.g [Tj 
[TTl not [36], multiplicative coalescent arises as the limit object in a large number of random graph 
models at criticality. Thus asymptotics of large random graphs has intimate connections with the 
general theory of coalescent (and fragmentation) processes that is currently a very active area in 
probability, see e.g. [5t [9l [38]. 

(d) Starting from an arbitrary configuration: In this work we have only considered the case where 
we start with the empty configuration 0„. One can imagine starting from a different configuration 
and then attempt to analyze the emergence of the giant component. Along the lines of [6], under the 
assumption that starting configuration satisfies suitable regularity conditions, it should be possible 
to study the asymptotic behavior of this model in terms of the entrance boundary of the standard 
multiplicative coalescent. This will be pursued in future work. 

(e) Proof techniques: The study of critical random graphs have generated an enormous amount 
of interest in the probabilistic combinatorics community and a large number of techniques have 
been developed for the fine scale asymptotics of such models, ranging from generating function and 
counting arguments (see e.g. [25], [E]); branching processes exploring local neighborhoods of the 
graph (see e.g [28j,and the references therein, [Mj); multi-type branching processes (see e.g. [T6] 
and the references therein); differential equation based techniques (see e.g. [H], [12] and for a 
comprehensive survey [32]); and the the breadth first walk approach coupled with the martingale 
central limit theorem (see e.g.|30j.[33] for one of the first few studies using this technique; [5] 
where this was used to explore the fine scale structure of the classical Erdos-Renyi random graph 
at criticality and connections with the multiplicative coalescent; see [36], [11] for further results 
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using this technique). Our methods are inspired by [7] which uses estimates on the size of largest 
component in the subcritical window to analyze the asymptotic behavior of the the sum of squares 
and cubes of component sizes near criticality. Estimates on the latter in turn allow the verification 
of the sufficient conditions given in [4j for convergence to the standard multiplicative coalescent. 
Although our general outline of the proof is similar to [7j as described above, it turns out that for 
general Achlioptas processes, completing this program is a rather challenging problem and for the 
BF model treated in the current work we need to develop quite a bit of machinery, which is done 
in Sections [5] and [U in order to obtain the required estimates. 

(f) Other structural properties of components: In this study we focus on the sizes of the compo- 
nents in the critical scaling window. Developing probabilistic methodology for understanding the 
actual structure of the large components is also of great interest, see for example [T7| for one of 
the first studies at a non-rigorous level, exploring the connection between the structure of these 
components at criticality and the internal structure of the minimal spanning tree in various random 
graph models (the "strong disorder" regime in statistical physics). In particular, it was conjectured 
in [17] that the diameter of the minimal spanning trees in various random graph models scales like 
n^/^. Rigorous studies have now been carried out for the Erdos-Renyi random graph model and 
fascinating connections have been discovered between the critical random graphs and the famous 
continuum random tree of Aldous (see [3], [2]). Proving such structural convergence of the entire 
components in dynamic network settings such as Achlioptas process at criticality, to random fractals 
such as continuum random trees, would be of great interest in a number of different fields. 



4 Proof idea 

Let us now give an idea of the proof. We begin by showing in Proposition 14. II below that, just before the 
critical window, the configuration of the components satisfies some important regularity properties. This 
proposition will be used in Section [8] in order to apply a result of [4] that gives sufficient conditions for 
convergence to the multiplicative coalescent. 

Proposition 4.1 Let 7 G (1/6, 1/5) and define tn = tc — n"'^ . Then we have 



^4/3 ^-7+1/3 



Now note that tn can be written as 



where 



52 (in) a 

n'/^CL'\tn) 

•S2{tn) 



tn = tc + ^ 



/? (4.1) 
(4.2) 
0. (4.3) 



n 



1/3 



„-7+l/3 

as n — )• 00. The above proposition implies that the configuration of rescaled component sizes, for large n 
at time "—00" , satisfy the regularity conditions for the standard multiplicative coalescent (see Proposition 
4 in [4J). 

Once the above has been proved, the second step is to show that through the critical window, the 
component sizes merge as in the multiplicative coalescent, at rate proportional to the product of the 
rescaled component sizes. This together with arguments similar to [7] will complete the proof of the main 
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result. 



Let us now outline the framework of the proof: 

• The bound on the largest component Cn ' (i) when t t (Proposition IS.ip plays a crucial role in 
proving the statements in Proposition 14.11 In order to achieve this, we introduce a series of related 
models from Section [5.21 through Section [5.4[ 

• Section[6]uses these models to prove asymptotically tight bounds on the size of the largest component 
through the subcritical window. The main goal of this Section is to prove Proposition 15. 1[ 

• Section [7] uses the bounds on the largest component from Proposition 15.11 to analyze the sum of 
squares and cubes of component sizes near the critical window. As one can imagine when working 
so close to the critical window, one needs rather careful estimates. Our arguments are based on 
a precise analysis of quadratic variation processes for certain martingales associated with the BF 
model. This section completes the proof of Proposition 14.11 

• Finally, in Section [8] we use Proposition 14. II and a coupling with the standard multiplicative coales- 
cent, in a manner similar to [7], to prove the main result. 

Without further ado let us now start with the proofs. 

5 An estimate on the largest component 

The following estimate on the largest component is the key ingredient in our analysis. Recall that tc 
denotes the critical time for the BF process. 

Proposition 5.1 Let 7 e (0,1/5) and let In{t) = Cn\t) be the largest component o/BF„(t). Then, for 
some B = ^(7) G (0, 00), 

P{/„(t) < m{n,t),yt < tc — n^'^} — )■ 1, when n — )■ 00 

where 

= b|^. (5.1) 

The proof of Proposition 15.11 will be completed in Section 16.31 In the current section we will give 
constructions of some auxiliary random graph processes that are key to our analysis. Although not 
pursued here, we believe that analogous constructions will be key ingredients in treatment of more general 
random graph models as well. The section is organized as follows. 

• In Section 15.11 we give the basic notation and mathematical conventions used in this paper. 

• In Section 15.21 we will carry out a preliminary analysis of the BF process and identify three deter- 
ministic maps ao,feo;Co from [0, cxd) to [0, 1] that play a fundamental role in our analysis. 

• Guided by these deterministic maps, in Section 15.31 we will define a random graph process with 
immigrating vertices and attachments (RGIVA) which is simpler to analyze than, and is suitably 
'close' to, the Bohman-Frieze process. A precise estimate on the approximation error introduced 
through this model is obtained in Section [6.31 
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• In Section 15.41 we will introduce an inhomogeneous random graph (IRG) model associated with a 
given RGIVA model such that the two have identical component volumes at all times. This allows 
for certain functional analytic techniques to be used in estimating the maximal component size. We 
will also make an additional approximation to the IRG model which will facilitate the analysis. 

• In Section 15.51 we summarize connections between the various models introduced above. 



5.1 Notation 

5.1.1 Graphs and random graphs 

A graph G = {V,£} consists of a vertex set V and an edge set £, where V is a subset of some type 
space X and is a subset of all possible edges {{^1,^2} : vi ^ V2 £ V}. An example of a type space is 
[n] = {1, 2, n}. Frequently we will assume X to have additional structure, for example to be a measure 
space {X,T,fi). When V is a finite set, we write |V| for its cardinality. 

G is called null graph if V = 0, and we write G = 0. G is called an empty graph if |V| = n and 
£" = 0, and we write G = 0„. 

Given two graphs, Gj = {Vi,£i} for i = 1,2, Gi is said to be a subgraph of G2 if and only if Vi C V2 
and £1 C £2 and we denote this as Gi < G2 (or equivalently G2 > Gi). We write Gi = G2 if Gi < G2 
and Gi > G2. 

A connected component C = {Vo,£o} of a graph G = {V,£} is a subgraph which is connected (i.e. 
there is a path between any two vertices in C). The number of vertices in C will be called the size of the 
component and frequently we will denote the size and the component by the same symbol. 

Let G be the set of all possible graphs {V,£) on a given type space X. When V is countable, we 
will consider Q to be endowed with the discrete topology and the corresponding Borel sigma field and 
refer to a random element of ^ as a random graph. All random graphs in this work are given on a fixed 
probability space which will usually be suppressed in our proofs. 

5.1.2 Probability and analysis 

All the unspecified limits are taken as n — )■ +00. Given a sequence of events {En}n>i, we say En (or E) 
occurs with high probability (whp) if ¥{En} — >• 1. For functions /, 5 : N — t- R, we write g = 0{f) if for 
some C £ (0,oo), lira snp g{n) / f (n) < C and g = ©(/) if 5 = 0(f) and / = 0{g). Given two sequences 
of random variables {^„} and {Cn}, we say ^„ = 0(Cn) whp if there is a C G (0, 00) such that ^„ < C(n 
whp, and write ^„ = &{Cn) whp if there exist < Ci < C2 < 00 such that CiCn < S,n < C'2Cn whp. 
Occasionally, when clear from the context, we suppress 'whp' in the statements. 

We also use the following little o notation: For a sequence of real numbers g{n), we write g = o{f) 
if limsup |g'(n)//(n)| = 0. For a sequence of random variables ^n, we write "^^ = Op{f)" if S,n/f{n) 
converges to in probability. 

For a real measurable function on a measure space {X,T,fJ-), the norms ||V'||2 and HV^Hoo are defined 

p d 

in the usual way. We use — > and — >• to denote the convergence in probability and in distribution 
respectively. 

We use =d to denote the equality of random elements in distribution. Suppose that {S,S) is a 
measurable space and we are given a partial ordering on S. Given two S valued random variables ^1,^21 
we say a pair of S valued random variables , ^2 given on a common probability space define a coupling 
of (^1,^2) if S,i =d £,* , i = 1,2. We say the S valued random variable ^1 stochastically dominates ^2, 
and write ^1 >d ^2 if there exists a coupling between the two random variables, say and such that 
Ci > C2 a.s. 
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For two sequences of S valued random elements ^„ and we say "^^ <cZ (,n whp." if there exist a 
coupling between ^„ and ^„ for each n (denote as ^* and S,*) such that ^* < ^* whp. 

Two examples of S that are relevant to this work are P([0,T] : M) and P([0,T] : Q) with the natural 
associated partial ordering. 



5.1.3 Other conventions 

We always use n,m,k,i, j to denote non-negative integers unless specified otherwise. We use s,t,T to 
denote the time parameter for continuous time (stochastic) processes. The scaling parameter is denoted 
by n. Throughout this work T = 2tc which is a convenient upper bound for the time parameters of 
interest. 

We use di,d2,-.- for constants whose specific value are not important. Some of them may appear 
several times and the values might not be the same. We use Ci,C2, ... for constants that appear in the 
statement of theorems. 



5.2 A preliminary analysis of Bohman- Frieze process 

Recall that BF„(t) denotes the BF process at time t and note that BF„ defines a stochastic process 
with sample paths in I?([0,T] : Q). Also recall that Cn\t) denotes the size of the i*'* largest component 
in BF„(t), Cn{t) = {Cn\t) : i > 1) is the vector of component sizes and Xn{t) denotes the number of 
singletons in BF„(t). We let J-t = T't = cr{BF„(s), s <i} and refer to it as the natural filtration for the 
BF process. 

At any fixed time t > 0, let COM{t) denote the collection of all non-singleton components 

CO>[(t) = {C«(t):|C«(t)|>2}. 

Recall that x{t) = Xn{t)/n. We will now do an informal calculation of the rate at which an edge 
e = {vi,V2} is added to the graph BF(t). There are three different ways an edge can be added: (i) both 
vi and V2 are singletons, (ii) only one of them is a singleton, (iii) neither of them is a singleton. 



Analysis of the three types of events: 



(i) Both vi and V2 are singletons. We will refer to such a component that is formed by connecting two 
singletons doubleton. This will happen at rate 

'Mn.al{x{t)). (5.2) 

The first product in the squared brackets is the count of all possible e = (61,62) S such that ei 
joins up two singletons and thus will be added to the graph, while the second product is the count of 
all e = (61,62) G £^ such that the first edge 61 does not connect two singletons while 62 connects two 
singletons and will be added. 
Define : [0, 1] [0, 1] as 

«ofe) = 2(f i+(i4)f)=i(.^ + (l-.V). (5.3) 

It is easy to check that 

On(^(*)) = CLQ{x{t)) +ra{t), where, sup|ra(t)| < 5/n. (5.4) 

t 

Recall that x{t) is the solution of the differential equation (12. 2p . To simplify notation we will write 
a*(x(t)) = a*{x) = a*{t) = a* and ao{t) = ao(x(t)) exchangeably. Similar conventions will be followed 



x„it) 

2 



+ 



Xn{t) 

2 



Xn{t) 

2 
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for the functions c* , cq and 5* , bo that will be introduced below. We shall later show that supj<7^ \xn{t) — 
x{t)\ — )• in probability (see Lemma [6.41 of this paper, also see t^lj)- This in particular implies that 
supj<j. |a* (t) — ao(t)| — ?• in probability. 

(ii) Only one of them is a singleton: This will happen if and only if ei does not connect two singletons 
while 62 connects a singleton and a non-singleton, thus at the rate 



(n-X,(t))X„(t). (5.5) 



We are also interested in the rate that a given non-singleton vertex (say, vq) is connected to any singleton 
which is 



3U9/ \ 9 . . X„(t) =^ <(x(t)). (5.6) 



Thus at time t a singleton will be added to COM{t) during the small time interval {t,t + dt], by attaching 
to a given vertex vq £ COM.{t), with the rate c*{t). 
Define cq : [0, 1] [0, 1] as 

co(y) = (i-y')y, 2/G [0,1]. (5.7) 

Then 

c*{x{t)) = Co{x{t)) +rc{t) and sup|rc(t)| < 2/n. (5.8) 



(iii) Neither of them is a singleton: This will happen at the rate 

J^^n(t)^^ fn-Xnit) 

Also, the event that two fixed non-singleton vertices are connected has the rate 



(5.9) 



3 M 9 / \ 9 / / - Kim)- (5.10) 
Let 6o : [0, 1] — ?• [0, 1] be defined as 

6o(y) = l-?/', yG[0,l]. (5.11) 

Then 

h*{x(t)) = bo{x{t)) + n{t) and sup \n{t)\ < 2/n. (5.12) 

t 

Note that for the study of the largest component one may restrict attention to the subgraph COA4{t). 
The evolution of this subgraph is described in terms of stochastic processes a* {x{t)),b* {x{t)) and c*{x{t)). 
In the next subsection, we will introduce a random graph process that is "close" to COM{t) but easier to 
analyze. Intuitively, we replace a*{t), b*{t),c*{t) with deterministic functions a{t), b{t),c{t) which are close 
to ao{t),bQ{t), co{t) (and thus, from Lemma [6^ whp close to a*{x{t)),b*{x(t)),c*{x{t))) and construct a 
random graph with similar dynamics as COM{t). 



5.3 The random graph process with immigrating vertices and attachment 

In this subsection, we introduce a random graph process with immigrating vertices and attachment 
(RGIVA) . This construction is inspired by [7] where a random graph with immigrating vertices (RGIV) is 
constructed - we generalize this construction by including attachments. RGIVA process will be governed 
by three continuous maps a, b, c from [0, T] — )• [0, 1] (referred to as rate functions) and the graph at time 
t will be denoted by IA„(t) = IA„,(a, 6, c)^ . When (a, 6, c) is sufficiently close to {ao,bo,co) , the RGIVA 
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model well approximates the BF model in a sense that will be made precise in Section 16.31 



The RGIVA process IA„(t) = IAn{a,b,c)t- Given the rate functions a,b,c, define IA„(t) as 
follows: 

(a) IA„(0) = 0, the nuh graph; 

(b) For t £ [0,T), conditioned on IA„(t), during the small time interval {t,t + dt], 

• (immigration) a doubleton (consisting of two vertices and a joining edge) will be born at rate n-a{t), 

• (attachment) for any given vertex vo in IA„(t), a new vertex will be created and connected to vq 



• (edge) for any given pair of vertices vi,V2 in IA„(f), an edge will be added between them at rate 



The events listed above occur independently of each other. 

In the special case where a{t) = b{t) = 1, c{t) = 0, and doubletons are replaced by singletons, the 
above model reduces to the RGIV model of [7]. We note that the above construction closely follows our 
analysis of three types of events in Section 15.21 replacing stochastic processes a* (x„ (t ) ) , (x„ (t ) ) , c* (x„ (t ) ) 
with deterministic maps a{t),b{t),c{t). 

The following lemma establishes a connection between the Bohman-Frieze process and the RGIVA 
process. Recall the partial order on the space P([0, T] : G). 



Lemma 5.2 Let {aL,bL,CL) and {au,bu,cu) be rate functions. Further, let U = Un be the event that 
{a*{t) < au{t),b*{t) < bu{t),c*{t) < cu{t) for all t e [0,r]} and L = Ln be the event that {a*{t) > 
aL{t),b*{t) > bL{t),c*{t) > CL{t) for all t e [0,T]}. Define forte [0, T] 



Then 

(i) Upper bound: COM^ <d lA^ = lAn{au,bu,cu). 

(ii) Lower bound: COM^ >d lA^ = IA„(ai, 5/,, c/,). 

Proof: We only argue the upper bound. The lower bound is proved similarly. Construct lA^(t) iter- 
atively on [0,T] as described in the definition, and construct COM^it) simultaneously by rejecting the 
proposed change on the graph with probabilities (1 — a*/af/)+, {1 — b* /bu)^ and (1 — c* /cu)'^ according 
to the three types of the events. Let r = inf{0 < t < T : a*{t) > au{t) or b*{t) > bu{t) or c*{t) > cu{t)} 
and set COM^{t) to be the null graph whenever t > t. This construction defines a coupling of lA^ and 
COMH such that COMH < IA?f a.s. The result follows. ■ 



5.4 An inhomogeneous random graph with a weight function 

In this section we introduce a inhomogeneous random graph (IRG) associated with IA„(a,6, c) for given 
rate functions a, b, c. For a general treatment of IRG models we refer the reader to [16], which our pre- 
sentation largely follows. We generalize the setting of [16] somewhat by including a weight function and 
considering the volume of a component instead of the number of vertices of a component. We begin with 
a description and some basic definitions for a general IRG model. 

A type space is a measure space (A',T, ^) where ^ is a complete separable metric space (i.e. a 
Polish space), T is the Borel u-field and ^ is a finite measure. 



at rate c{t) 




on [/' 
on U 




on 
on L 
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A kernel on the type space [X^T^ii) is a measurable function k : X x X ^ [0, oo). The kernel k is 
said to be symmetric if k(x, y) = ^(y, x) for all x, y G . We will also use x, y instead of x, y for elements 
in X when there is no confusion between an. x ^ X and the function x{t) defined in (12. 2p . 

A weight function (/> is a measurable, non-negative function on (Af,T, /u). 

A basic structure is a triplet {{X, T, fj.), k, (j)}, which consists of a type space, a kernel and a weight 
function. 

The IRG model: Given a type space {X,T,fi), symmetric kernels {nn}n>i, and a weight function 
(f>, a random graph R,Gn(Kn) ( = R-Gn(Kru/-i) = -E^Gn(Kra, /-t, <j))), for any integer n > 0, is constructed as 
follows: 

(a) The vertex set V are the points of a Poisson point process on {X, T) with intensity n ■ //. 

(b) Given V, for any two vertices x, y G V, place an edge between them with probability • Kn{x, y)) A 1. 

One can similarly define an IRG model associated with a basic structure {{X jT, n), K,(p}, where k is 
a symmetric kernel, by letting Kn = k for all n in the above definition. 

The weight function (j) is used in defining the volume of a connected component in the above construc- 
tion of a random graph. Given a component of RG„(k, ^u, (/>) whose vertex set is Vo, define J2xeVo '?^(^) 
as the volume of the component. 

One can associate k with an integral opertor /C : L'^ip) — )• L'^ip) defined as 

/C/(x)= / K{x,y)f{y)ii{dy) (5.13) 
Denote by p = p(k) the operator norm of /C. Then p = p{i<i) = ||/C|| = supujn^^i ||/C/||2- 

Given rate functions a, 6, c, there is a natural basic structure and the corresponding IRG model asso- 
ciated with IA„(a, 6, c), which we now describe. 

Fix t G [0,T]. Then the following two stage construction describes an equivalent (in law) procedure 
for obtaining IA„(a, 6, c)t ■ 

Stage I: Recall that transitions in IA„(a,6, c) are caused by three types of events: immigration, 
attachment (to an existing vertex) and edge formation (between existing vertices). Consider the random 
graph obtained by including all the immigration and attachment events until time t but ignoring the 
edge formation events. We call the components resulting from this construction as clusters. Note that 
each cluster consists of exactly one doubleton (which starts the formation of the cluster) and possibly 
other vertices obtained through later attachments. Note that doubletons immigrate at rate a(s) and 
supposing that a doubleton is born at time s, the size of the cluster at time s < u < t denoted by w{u) 
evolves according to a integer-valued time-inhomogeneous jump Markov process starting at w{s) = 2 and 
infinitesimal generator A{u) given as 

A{u)f{r) = c{u)r ■ {f{r + 1) - /(r)) , / :N^M,s<n<t. (5.14) 

We set w{u) = for < n < s and denote this cluster which starts at instant s by {s,w). 

Stage II: Given a realization of the random graph of Stage I, we add edges to the graph. Each pair of 
vertices will be connected during (s, s+ds] with rate ^h{s). Thus the number of edges between two clusters 
x = {s,w),y = {r,w) at time instant i is a Poisson random variable with mean ^ w{u)w{u)b{u)du. 
Consequently, 

1 /■* 

P{x and y is connected | Stage 1} = 1 — exp{ / w{u)w{u)h{u)du} (5.15) 

n Jo 

1 /■* 

< — / w{u)w{u)h{u)du. (5.16) 
n Jo 
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It is easy to see that the graph resulting from this two stage construction has the same distribution 
as IAn{a,b,c)t. 

We now introduce an IRG model associated with the above construction in which each cluster is treated 
as a single point in a suitable type space and the size of the cluster is recorded using an appropriate weight 
function. Let X = [0, T] x W, where W = P([0, T] : N) is the Skorohod -D-space with the usual Skorohod 
topology. Denote by T the Borel sigma field on [0, T] x W. For future use, we will refer to this particular 
choice of type space {X, T) as the cluster space. For a fixed time t > 0, consider a weight function defined 
as 

0t(x) = u;(t), x= G [0,r] X W. (5.17) 

Then this weight function associates with each 'cluster' x its size at time t. We now describe the finite 
measure that governs the intensity of the Poisson point process 'Pt(a, b, c) of clusters (regarded as points 
in X). Denote by Ug the unique probability measure on the space W under which, a.s., w{u) = for 
all u < s, w{s) = 2 and w{u),u G [s,T] has the probability law of the time inhomogeneous Markov 
process with generator {A{u),s < u < T} defined in (I5.14p . Let fi he a finite measure on X defined as 
li{dsdw) = i's{dw)a{s)ds, namely, for a non-negative real measurable function f on X 

f{x)dfj.{x) = / a(s) I / f{s,w)dvs{w) I ds. 



X Jo \Jw / 

We also define for each t G [0,T], a finite measure fit on X by the relation fit{A) = /"(^ H ([0, t] x W)). 
Then for / as above, 

/ /(x)(i/it(x) = / a{s)( [ f{s,w)dus{w)]ds. (5.18) 
Jx Jo \Jw / 

The measure fit will be the intensity of the Poisson point process on X which will be used in our con- 
struction of the IRG model associated with IA„(a, b, c)t. Now we describe the kernel that will govern the 
edge formation amongst the points. Define 



Kn,t(x, y) = Kn,t((s, w), (r, w)) = n \ 1 - exp{ / w{u)w{u)b{u)du] . (5.19) 



n 







We will also use the following modification of the kernel Kn^t- 



Kt{-x,y) = Kt{{s,w),{r,w)) = / w{u)w{u)b{u)du. (5.20) 

Jo 

With the above definitions we can now define IRG models RG„(Kn,t, fJ-t-,<Pt) and RG„(Kt, ^t, (pt) associated 
with the type space {X,T,fJ-t)- 

Denote the size of the largest component [resp. the component containing the first immigrating 
doubleton] in IA„(a,6, c)t by C^^\a,b,c)t [resp. C^'^\a,b,c)t]. Also, denote the volume of the largest 
component [resp. the component containing the first cluster] in HGniKt, fJ,t,(pt) by C^^^Kt, iJ,t,(pt) [resp. 
C^^") {Kt, i^t, (t>t)]- Define C(^)(Kn,t, A^t, </>t), C^^\Kn,t, l^tj^Pt) in a similar fashion. The following is an immedi- 
ate consequence of the above construction. 

Lemma 5.3 We have 

(C(i) (a, 6, c)t, C(°) (a, b, c)t) =d (C^^) {Kn,t,fit, (t>t),C^''^ {'^n,t,fit, (t>t)) 

and 

C'-^\Kn,t,fJ-t,(t>t) <dC''^HKt,fJ-t,4>t), C'-°\Kn,t,fJ-t,(t>t) <dC^^\Kt,fJ.t,(l)t)- 

For future use we will write RG„(kj, ^t, cpt) = RG„^j(a, 6, c). 
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5.5 A summary of the models 



As noted earlier, the key step in the proof of Proposition 14.11 is a good estimate on the size of the largest 
component in the Bohman-Frieze process BF„(t) as in Proposition 15. 1[ For this we have introduced a 
series of approximating models. We summarize the relationship between these models below. 

• We can decompose the Bohman-Frieze process as BF„ = COA4n U namely the non-singleton 
components and singleton components at any time t. 

• We shall show that COM.n ~ IA„(ao, bo, cq), where oq, bo, cq are defined in (j5.3p . (|5.1ip . ()5.7p . More 
precisely we shall show that as n — )• oo, for any fixed 6 > 0, we have, whp. 

IA„((ao - 6)+, {bo - 6)+, (co - 6)+) <d COMn <d IA„((ao + 6)A1, {bo + 6)Al, (co + 6) A 1). 

This is a consequence of Lemma 15. 2[ 

• Given rate functions {a,b,c), for all t G [0,T], 

C«(a,6,c)t =rfC«(K„,t,/it,0O <dC(')(/€t,/xt,0t), i = 0,l. 

Here Kn,t, i^t, fJ-tjfpt and a,b,c are related through (I5.19|) . ()5.20p . (I5.18P (see also ()5.14ll ). ()5.17p . 
respectively. 

6 Analysis of the largest component at sub-criticality: Proof of Propo- 
sition 15.11 

This section proves Proposition 15.11 The section is organized as follows: 

• In Section [6.11 we reduce the problem to proving Proposition 16.31 We give the proof of Proposition 
15.11 using this result. Rest of Section [6] is devoted to the proof of Proposition 16.31 

• In preparation for this proof, in Section 16.21 we present some key lemmas that allow us to estimate 
the errors between various models summarized in Section 15.51 Proofs of some lemmas (Lemmas 16. 6^ 
16.91 and l6.10p are deferred to later sections. 

• Using these lemmas, in Section lOl we prove the key proposition. Proposition 16.31 The rest of Sec- 
tion [6] proves the supporting Lemmas 16.61 16.91 and 16.101 

• In Section 16.41 we introduce a branching process related to the IRG model, and prove Lemma 16.61 
A key step in the proof is Lemma 16.131 whose proof is left to Section 16.51 

• Section [6.61 analvzes the kernel Kt associated with RG,i^f (a, 6, c) and proves Lemma 16.91 . 

• Finally, in Section 16.71 we give the proof of Lemma 16.101 
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6.1 Prom the largest component to the first component 



In this section we will reduce the problem of proving the estimate on the largest component in Proposition 
15.11 to an estimate on the first component as in Proposition 16.31 This reduction, although somewhat 
different, is inspired by a similar idea used in [7]. 

Recall that Cn\t) = In{t) denotes the largest component in BF„(f). Let C^(t), < s < t, denote the 
component whose first doubleton is born at time s in BF„(t). In particular C^(t) = if there is no 
doubleton born at time s. Without loss of generality, we assume that the first doubleton is born at time 
0. Then C^{t) denotes the component of the first doubleton at time t of the BF process. The following 
lemma estimates the size of the largest component In{t) in terms of the size of the first component. 

Lemma 6.1 For any n G N, to G [0)^] c-i^d deterministic function a : [0,T] — t- [0,oo) 

P{/„(t) > a{t), for some t < to} < nrP{C°(t) > a{t), for some t < to}. 

Proof: Let {'BFn\t),t > OjjgNo be an i.i.d. family of {BF„(t),t > 0} processes on the same vertex set 
[n]. Let be a rate n Poisson process independent of the above collection. Denote by {Tj,i G N} the 

jump times of the Poisson process. Set tq = 0. Denote the first component of BF^*^ at time t by Jn\t). 
Consider the random graph 

and let /^(t) denote the size of the largest component in G^. Then since a* (t) < 1 for all t, /„ <ci 
Thus 

P{/„(t) > a{t), for some t < to} < P{/^(t) > a{t), for some t < to} 

= ^ P{/^(t) > a{t), for some t < to, N{T) = k} 
fceNo 

F{J^'\t) > a{t), for some t < to, for some i < k}F{N{T) = k} 

fcGNo 

kF{C^{t) > a{t), for some t < to}F{N{T) = k}. 

fceNo 

The result follows. ■ 



< 



< 



Next, in the following lemma, we reduce an estimate on the probability of the event {C^(t) > 
a{t), for some t < to} to an estimate on sup^gjo.to] '^i^)^{^ni^) > a{t)}. 

Lemma 6.2 There exists an No & N such that for all n > A'^o, to £ [0, T] and continuous a : [0, T] — )• 
[0,oo) 

P{C°(t) > 2a{t), for some < t < to} < IGnT^ sup {a(s)P{C°(s) > a{s)}} . (6.1) 

0<s<to 

Proof: Fix No £ N such that for all n > No, supg^^o,T]{o-ni^) V b^is)} < 2. Consider now n > Nq. Define 
T = inf{t > : C^{t) > 2a{t)}. Then 

P{C°(t) > 2a{t) for some t G [0, to]} = P{r < to}. (6.2) 

Denote by -H-t the event that components and merge at time t. By convention this event is 
taken to be an empty set if no doubleton is born at time instant s. Then 

{T = t} = {CO(t_) < 2a(t)} n {CO(t-) +C^(t-) > 2a(t);C0 for some s < t}. 

Next note that 
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• Since a* (s) < 2, the rate at which doubletons are born can be bounded by 2n. 

• Given a doubleton was born at instant s, the event {C^ C^, for some u G {t,t + dt]} occurs, con- 
ditionahy on Tt, with probabihty ^C^^{t)C:^{t)b^{t)dt. This probabihty, using the fact that b^{s) < 2 
and C*(t) < n, on the event {C^{t) < 2a{t)} is bounded by 4a(t)di. 

• P{C°(t) + C^{t) > 2a{t)} is bounded by 2P{C°(t) >a{t)}. 
Using these observations we have the fohowing estimate 



[r<to} < E / l{co(t)<2a(t)} I na:(s) • (ic°(t)C:(t)) 
J[o,to] [J[o,t] n 

I / 2n • 2P(C°(t) > a{t)) ■ {Aa{t))ds 



dt 



< 



< 



[OM lJ[o,t] 

{2nt) ■ 2P(C°(t) > a{t)) ■ {Aa{t))dt 
[0,to] 



< IGnT^ sup {Q(t)P{C°(t) > a{t)]] . 

iG[0,to] 

Result fohows on combining this estimate with (j6.2p . ■ 

The following proposition will be proved in Section 16.31 

Proposition 6.3 Given rj £ (0, oo) and 7 G (0,1/5), there exist B,C,Ni G (0, 00) such that for all 
n>Ni 

P{C°(t) >m{n,t)/2} <Cn-'^ for alio <t <tc-n-^, (6.3) 
where m{n,t) is as defined in (jS.ip . 

Remark: Intuitively, one has that in the subcritical regime, i.e. when t < tc, P{C°(t) > m} < die-'^^m 
for some constants di,d2- This suggests a bound as in (|6.3p for each fixed t < tc- However, the constants 
di and d2 depend on t, and in fact one expects that, d2{t) — t- when t '[ tc- On the other hand, in order 
to prove the above proposition one requires estimates that are uniform for all t < tc — n~'^ as n — t- 00. 
This analysis is substantially more delicate as will be seen in subsequent sections. 



We now prove Proposition 15.11 using the above results. 



Proof of Proposition I5.lt Fix 7 G (0, 1/5) and fix r/ > 2 + 27. Let B, C, Ni be as determined in 
Proposition 16.31 for this choice of 7^,7 and let m{n,t) be as defined in (|5.ip . Without loss of generality 
we can assume that Ni > Nq where A'^q is as in Lemma 16.21 Then applying Lemmas 16.11 and 16.21 with 

to = tc — n~"' and a{t) = m{n, t), we have 

F{In{t) > m{n, t), for some < t < tc - n"^} < nT¥{C^{t) > m{n, t), for some < t < tc - n"^} 

< 16n^r^ sup {m{n,s)F{C^{s) >m{n.s)/2}} 

sG[0,tc-n-T] 

< 16C5n2-^+2Tr3(logn)^ 

Since ?? > 2 + 27, the above probability converges to as n — )• 00. The result follows. ■ 
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6.2 Some Preparatory Results 

This section collects some results that are helpful in estimating the errors between various models de- 
scribed in Section 15.51 

The first lemma estimates the error between Xn(t) = x{t) = Xn{t)/n and its deterministic limit x{t) 
defined in (lO). 



Lemma 6.4 For any T > 0, there exists a C{T) G (0, oo) such that, for all 71 G [0, 1/2), 

F{ sup \xn{t) - x{t)\ >-^}< exp{-C7(T)ni-2^^}. 

Proof: Recall that [n] = {1, 2, . . . Let En = n~^[n] and let E = [0, 1]. Recall the three types of 
events described in Section 15.21 that lead to edge formation in the BF model. Of these only events of type 
(i) and (ii) lead to a change in the number of singletons. For the events of type (i), i.e. in the case when 
a doubleton is created, x decreases by 2/n. Two key functions (see (|5.2p ) for this case are 

/-2(y) = <iy) 

/_2(y) = ao(y) = ^(y' + (i-y')y). 

For the events of type (ii), i.e. in the case when a singleton attaches to a non-singleton component, x 
decreases by 1/n. Two key functions (see (j5.5p l for this case are 

/Ii(y) = (l-yK(y) 

/_i(y) = (1 - y)co(y) = y(l - y^)il - y). 

Note that < fi{x) < 1 for / = — 1, —2, and that x{t) is a Markov process on the state space En for 
which at time t we have the transitions x{t) ~^ x{t) — 1/n at rate nfli{x{t)) and x{t) — ^ x{t) — 2/n at 
rate n/*2(a;(t)). Furthermore 

l/Ii(y) - f-iiy)\ < - \fl2iy) - /-2(y)l < -, for all y e [O, l]. (6.4) 
n n 

Let y_i(-), y_2(-) be independent rate one Poisson processes. Then the process x{t) started with a;(0) = 1 
can be constructed (see eg. [32] , ) as the unique solution of the stochastic equation 

1„ / /■* \ 2_ ^ 



x{t) = 1 - -y_i r_^ixis))dsj - -y_2 fl^{x{s))ds j . (6.5) 

By Equation (j2.2p . the limiting function x(-) is the unique solution of the integral equation 

x{t) = 1- f f-i{x{s))ds - f 2f^2{x{s))ds. (6.6) 
Jo Jo 

Also note that ^y, z £ E 

|(/_i(y) + 2/_2(y)) - {f-i{z) + 2/_2(^))| < Q\y - z\. (6.7) 
Using (16. 6p and (16. 5p we get 

\xit)-xit)\<A^,{t) + A^{t)+A^{t) 
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where 



Am 

and by ([63]) 



l=-l-2 



-Yi (n f fl{x{s))ds 
n \ Jo 



< 4 sup sup 

l=-l-2 t<T 



Yiint) 



n 



Am 



^ l[fnx{s))-fiix{s))]ds 



Z=-l,-2 



<It. 

n 



and finally by ([67 



Am 



Yl nMxis))-Mx{s))]ds 



Z=-l,-2 



< 6 / — x(s)|(is. 

'o 



Combining these estimates we get 



/ 7 

\x{t) — x{t)\ < — 1-4 sup sup 



Yiint) 



n 



+ 6 / — x{s)\ds. 

/o 



This implies, by Gronwall's lemma (see e.g. [23], p498) 
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sup |x(s) — x(s)| < — 1-4 sup sup 

s<T \n /=_! 2t<T 



Yiint) 



n 



o6T 



Proof is completed using standard large deviations estimates for Poisson processes. ■ 

In the next lemma we note some basic properties of the integral operator associated with a kernel k 
on a finite measure space. 

Lemma 6.5 Let k, k' be kernels given on a finite measure space (X, T, /i). Assume that k, k' £ L'^il-i x fi). 
Denote the associated integral operators by K, and K! (see (I5.13P ) and there norms by p{k),p{k') respec- 
tively. Then 

1/2 



(i) K, is a compact operator. In particular p{k) = ||/C|| < ||k||2 
(a) If K < k' , then p{k) < p{k'). 

(Hi) p{k + k') < p{k) + p{i^') and p{tK) = tp{K) for t > 0. 

(iv) \p{k)-p{k')\<p{\k-k'\). 

(v) p{k) < \\k\\ooP{X)- 



1^ {x,y)p{dx)p{dy) 



< oo. 



Proof: (i) is a standard result, see Theorem VI. 23 of 
(ii) For any nonnegative / in L^{p), K,f{x) < K,' f{x) pointwise. Thus for such /, ||/C/||2 < ||/C'/||2. 
Result follows on observing that the suprema of ||/C/||2, ||/C'/||2 over {/ € : ||/||2 = 1} is the same as 
the suprema over {/ G : ||/||2 = 1, / > 0} . 

(ifi) This follows immediately from the facts that ||(/C + /C')/||2 < W^fh + W fh and /C(t/) = tKf. 
(iv) Note that k < k! + \k — k'\ and n' < k + \k — Result follows on combining this observation with 
(ii) and (iii). 
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(v) This follows immediately from (i) and the fact that ||k||2 < ||«^||ooM(^Y). 



We now present some auxiliary estimates for the IRG model from Section [5.41 The following lemma 
will be proved in Section 16.41 Recall the definition of a basic structure from Section 15.41 

Lemma 6.6 Let {(A', 7", /i), k, </>} be a basic structure, where k is symmetric. Suppose that fj, is non- 
atomic and p{k) = ||/C|| < 1. For fixed xq G X, denote by C^^{xo) the volume of the component of 
RG„,(k) that contains xq. Define C^'^(xo) = if xq is not a vertex in RG„(k). Then for all m £ N 

P{C^^(xo) > m} < 2exp{-CiA2m} (6.8) 

where 



The above result will be useful for estimating the size of a given component in RG„(a, b, c). One difficulty 
in directly using this result is that the kernel Kt and the weight function cpt defined in (15.20p and (|5.17p 
are not bounded. We will overcome this by using a truncation argument. In order to control the error 
caused by truncation, the following two results will be useful. For rest of this subsection the type space 
[X^T) will be taken to be the cluster space introduced above (j5.17p . 

Lemma 6.7 Given rate functions {a,b,c) and t £ [0,T], let fit be the finite measure on {X,T) defined 
as in (jS.lSp . Let Vn be a Poisson point process on {X,T) with intensity n ■ fj,f Define 

Yn sup wit). 

Then for every A £ (0, cxd) 

F{Yn > A} <2T- n{l - e"^)^/^. 

Proof: Let N be the number of points in "P^, then N is Poisson with mean f^na{s)ds < uT. Let 

{Z2*''}i>i be independent copies of Z2 (also independent of A^), where Z2 is a pure jump Markov process 
on N with initial condition ^2(0) = 2 and infinitesimal generator ^0 defined as 

A^fik) = k{f{k + 1) - f{k)), ken, f-.n^m. 

Thus Z2 is just a Yule process started with two individuals at time zero. Note that 

Yn< sup w{T) <d sup z!^'\t), 

{s,w)eV„ l<i<N 

where the first inequality holds a.s and the second inequality uses the fact that c < 1. Standard facts about 
the Yule process (see e.g. [37]) imply that Z2\t) is distributed as sum of two independent Geom{e~"^}. 
Thus 

P{y„ > ^} < E(iV) • ¥{Z2{T) > A} 
<nT-2{l-e~^)^/^. 

This completes the proof of the lemma. ■ 

The following corollary follows on taking A = C log n in the above lemma. 
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Corollary 6.8 LetYn be as in the above lemma and fix t] G (0,oo). Then there exist Ci{r]),C2{i]) € (0, oo) 
such that for any rate functions (o, b, c) 

¥{Yn > Ci(r?) logn} < C2{v)n"'^, for all n G N. 

From Section [5.21 recall the definitions of the functions ao,bo,CQ associated with the BF model. The 
following lemma will allow us to argue that RG„(ao, boi cq) is well approximated by RG„(a, 6, c) if the 
rate functions (a, 6, c) are sufficiently close to {ao,bo,co). Let {{X ,T, Ht), ^t, (j)t) be the basic structure 
associated with rate functions (a,6, c). Let ICt be the integral operator defined by ()5.13p . replacing 
there by {fit, Kt). Let pt = p{Kt). In order to emphasize the dependance on rate functions (o, b, c), we will 
sometimes write pt = pt{a,b,c). Similar notation will be used for Kt,pt,(pt and ICt- 

Lemma 6.9 Fix rate functions {a,b,c). Suppose that infgg[o,T] ciis) > and for some 6 G (0,oo), c(s) > 
Os, for all s G [0,r]. Given 6 > and t € [0,T], let 

p+,t = Pt{{a + <5) A 1, (6 + J) A 1, (c + J) A 1), p^^t = Pt{{a - 6)+, (b - 6)+, (c - 6)+). 

Then there exists C2 G (0, 00) and 5q G (0, 1) such that for all 5 < 5q and t G [0, T] 

Tn?ai{\pt- p+,t\,\pt- P~,t\] < C2(-log(5)^(5^/^ 

The proof of the above lemma is quite technical and deferred to Section 16.61 

The next lemma gives some basic properties of pt{ao,bo,co). Recall that tc denotes the critical time 
for the emergence of the giant component in the BF model. 

Lemma 6.10 Let p{t) = pt(ao, 69, cq). Then: 
(i) p{t) is strictly increasing in t £ [0,T]; 
(a) p{tc) = 1; 

(Hi) limg^o+(p(*c) - Pitc - s))/s = p'_{tc) > 0. 
The proof of the lemma is given in Section 16. 7i 

6.3 Proof of Proposition 16.31 

This section is devoted to the proof of Proposition 16. 3i Fix ry G (0,oo) and 7 G (0, 1/5). 
Step 1: from BF„ to IA„ 5 

Let 7i = 2/5 and define = {supo<t<T — x{t)\ < n"'^^}. 

From Lemma 16.41 

P{K^} < exp{-C(r)ni~27i} = exp{-C7(T)ni/5}. (6.10) 
Prom ()5.4p and recalling that the Lipschitz norm of oq is bounded by 2 (see ()5.3p ). we have that on En 

\a*{t) - ao{t)\ < 5n-^ + 2n~^\ for all t G [0,T]. 

Similar bounds can be shown to hold for b* and c*. Thus we can find ni G N and d\ G (0, 00) such that, 
for n> ni, on En 

a*n{t) < ao{t) + 6n,b*n{t) < bo{t) + 6n,c*n{t) < co(t) + <5„, for all t G [0,r], 



22 



where 5n = din '''i. Since a^,b*^,c^ are all bounded by 1, setting (ao(0 + Sn) A 1 = 0^,5 and similarly 
defining bn,s, Cn,s, we in fact have that 

<(0 < an4t),b*^{t) < bn,s{t),c*^{t) < Cn,5{t), for all t G [0,T]. 

Let C'/s(t) denote the size of the first component in IAji((2n,5i ^n,(5; Cn,(5)t. From Lemma 15.21 we have for 
any m G N 

P{CO(t) > m, En} < nCL%t) > m, E^} < IP{Q(t) > m}. (6.11) 

Step 2: from IA„ 5 to RG„,5,a 
For t G [0,T], and rate functions an^s,bn,s,Cn^s, consider the IRG model HGniKt^s^ (^-1,5, where Kt^s = 
i^t{(in,s-,bn,&-,Cn,s) and ^Jit,s is the measure for the IRG model corresponding to these rate functions as 
defined in ()5.18p . Let An = Ci{r])\ogn, where Ci{rf) is as in Corollarv 16.81 Consider the following 
truncation of the kernel Kt^s and weight function (j)t{s,w) = w{t): 

Kt,M(x,y) = Kt,5(x,y)l{^(r)<A„}l{^(T)<A„}, X = {s,w),y = {r,w) 

and 

(j)t,Ais,w) = <j)tis,w)l{^(^T)<A„}- 

Then ||(/'i,A||oo < <TAl. 

Recall the Poisson point process 'Pt{a,b,c) associated with rate functions {a,b,c), introduced below 
(!5.17p and write Vt,5 = 'Ptian,5 , bn,5 , Cn^s) ■ Let Yn^s = sup(^ ,„)gp^^ w{T). From Corollary [O] 

F{Yn,5 > An} < C2{r])n~\ (6.12) 

Let C^Kt) = C^^{an,s,bn,s,Cn,s) be the volume of the 'first' component in RG„,t(an,<5, ^n,5, c„,5) = 
ILGnint^S, fj.t,6,<Pt)- Then from Lemma [5^ 

P{Q(t) >m}< nC^^sit) > m}. (6.13) 

Letting C^^ y^(t) denote the volume of the first component in RGn(Kt^5^A) Mt,(5i '^t,A)i namely the random 
graph formed using the truncated kernel. Then 

P{C^J(t) > m} < F{Yn,5 > An} + nCn${t) > m, Yn,5 < An} 

= P{K„,5 > An} + nCn%A{t) > m, Yn,S < An} 

< F{Yn,s>An} + nC^%At)>m} (6.14) 

Step 3: Estimating C^^^ 
We will apply Lemma \6M replacing {{X, T, /i), k, cp} there by {{X, T, fJ't,5), i^t,5,A-, (t>t,A}-, where t G (0, — 
n~^). From ()6.8p we have 

P{Cn,M(t) >m}< 2exp{-C7iA2m}, (6.15) 

where 



8||</>i,yl||oo(l + ^^^t,5,A\\oc^J't,5iX)) ' 

and A = 1 — p{nt,s,A)- We now estimate p{Kt,6,A)- Since K.t^s,A < by (ii) of Lemma [6.5) we have 
p{i^t,5,A) < p{i^t,5)- Note that rate functions (ao,6o,co) satisfy conditions of Lemma 16.91 Thus, recalling 
that 5n = din^^/^, we have from this result, that for some d2 G (0,oo), p{K,t,8) < p{^t) + (i2(log n)'^n~^/^, 
for all t <T. Here Kt = Kt{ao, bo, cq). 

Next, by Lemma [6.10|. there exists ^3 G (0,oo) such that p{Kt) < 1 — d^^tc — t) for all t G [0, ic)- 
Combining these estimates, we have for t < tc — n~'^, 



Pii<-t,5,A) <l-d3{tc-t) + d2(logn)V 



-1/5 
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Recalling that 7 E (0, 1/5) we have that, for some 712 S (ni,oo) and £ (0,oo), 

p{i^t,&,A) ^ 1 — d/i{tc — t), for all t € (0, tc — n~'^) and n > 722- 

Using this estimate in ()6.15p and recalling that ||(/>t,A||oo < H/^t.^./illoo < > have that for some 
rfs G (0,00) 

P{C5 .(t) > m} < 2exp{-— -^^^^(tc - tfm}, for ah m e N,t G (0,tc - n"'^) and n > n2. (6.16) 
' ' (log ny 

Step 4: Collecting estimates: 
Combining (fOOjl . (l6TT]l . (l6J3]l . (l6J2]l . (l6H]l and (fOBD . we have 

P{C°(t) > m} < P{S^} + P{y„,5 > + P{C^,^(t) > m} 

< e-C(r)nV5 ^ + 2exp{-d5^^^— ^m}. (6.17) 

(log nj'^ 



Finally, result follows on replacing m in the above display with 



rjQognf 



The following lemma will be used in the proof of Lemma I6.10[ We will use notation and arguments 
similar to that in the proof of Proposition 16.31 above. 

Lemma 6.11 Let (a,b,c) be rate functions. Fix t E [0,T]. Let Il^{t) denote the largest component in 
IAn{a,b,c)t. Suppose that pt{a,b,c) < 1. Then for some Cq £ (0, 00) 

P{/^-^(t) > Co(logn)^} ^ when n^oo. 

Proof: Let Cl^{t) be the first component of IA„(a, 6, c)j. Then an elementary argument (cf. proof of 
Lemma l6JJ shows that for m > 



P{Z^^(t) >m}< Tn¥{C'^^{t) > m}. 

By an argument as in (I6.14p . we have 

P{C(t) > m} < P{Cr (t) >m}< P{y„ > An} + nC^^Ait) > m}, 

where C^*^, Yn and C^*^ correspond to C^^, Yn^s and C^l introduced above in the proof of Proposition 
16.31 with {an^&,hn^s-,Cn,s) replaced with {a,b,c). From Corollarv 16.81 we can find di £ (0, 00) such that 
P(^n > dilogn) = 0(n~^). Let An = dilogn. Then, recalling that pt{a,b,c) < 1, we gave by Lemma IBTHl 
that, for some ^2 G (0,cx)), 

P{Cn,A(i) >m} < 2exp{-d2m/(logn)3}. 

Taking m = ^(logn)^, we have F{C^^^{t) > m} = 0(n~^). Combining the above estimates we have 
P{/^-^(t) > ^(logn)''} = 0(n~i). The result follows. ■ 



6.4 Proof of Lemma 16. 6t A branching process construction 

The key idea in the proof of Lemma 16.61 is the coupling of the breadth first exploration of components in 
the IRG model with a certain continuous type branching process. This coupling will reduce the problem 
of establishing the estimate in Lemma [6.6l to a similar bound on the total volume of the branching process 
(Lemma 16. 13p . We refer the reader to [16] where a similar coupling in a setting where the type space X 
is finite using a finite-type branching process is constructed. In this subsection we will give the proof of 
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Lemma 16.61 using Lemma I6.13[ Proof of the latter result is given in Section 16.51 



Throughout this section we will fix a basic structure {{X,T, fj,), n,4>}, where «; is a symmetric ker- 
nel, and a xq £ X. Let RG„(k) be the IRG constructed using this structure as in Section [531 We 
now describe a branching process associated with the above basic structure. The process starts in the 
0-th generation with a single vertex of type xq G X and in the A:-th generation, a vertex x will have off- 
spring, independently of the remaining k-th generation vertices, according to a Poisson point process on X 
with intensity k{x, y)^{dy) to form the [k + lY^ generation. We denote this branching process as BP(2;o). 

Denote by {^l'^^}^-^ C X the fc*^ generation of the branching process. Define the volume of the k-th 
generation as Gk = Ylu^i The total volume of BP(xo) is defined as G = G{xo) = YlT=o ^k- 

The following lemma, proved at the end of the section, shows that C^^{xq) is stochastically dominated 
by G(xo). 

Lemma 6.12 For all mi S N, 

nCn^ixo) > mi} < F{G{xo) > mi}. 

Next lemma, proved in Section [6.51 shows that the estimate in Lemma 16.61 holds with C^'^(xo) replaced 
by G{xo). 

Lemma 6.13 Suppose that p{n) = ||/C|| < 1. Then for all m gN 

F{G > m} <2exp{-CiA'^m} (6.18) 

where A and Ci are as in (16. 9p . 

Using the above lemmas we can now complete the proof of Lemma [6.61 

Proof of Lemma 16. 6t Proof is immediate from Lemmas 16.131 and 16.121 . ■ 

We conclude this section with the proof of Lemma 16.121 

Proof of Lemma 16.12t Without loss of generality assume that C^^{xo) ^ 0. We now explore the 
component C^'^{xo) in the standard breadth first manner. 

Define the sequence of unexplored sets {Um}m>o and the set of removed vertices {Rm}m>o 
iteratively as follows: Let Rq = 0, C/q = {^^o} and yi = xq. Suppose we have defined Rj,Uj, j = 
0, 1, • • • ,m-l and Um-i = {ym,ym+i, ■ ■ ■ Vt^}- Then set 

Rm = Rm-1 U {ym} 

Um = Um-1 U Em \ {ym} 

where 

Em = {x £ X : x is a neighbor of ym in RG„(k) (i.e. {x, ym} is an edge) and x ^ Rm-i U Um-i}- 

If Um-i = we set Uj = Ej = % and Rj = Rm-i for all j > m. Thus Um are the vertices at step m 
that have been revealed by the exploration but whose neighbors have not been explored yet. Note that 
the number of vertices in Rm-i U Um-i equals tm- Label the vertices in Em as yt™.+2, • • • ytm+\Em\- 

With this labeling we have a well defined specification of the sequence {Rj,Uj, Ej^i}j,=^g. Note that 
Q'^ixo) = "T-o if and only if Umo-i 7^ 0, Umo = and \Rmo \ = "t-o- 
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We will now argue that for every m S N, conditioned on {Um~i, Rm-i}, Em is a Poisson point process 
on the space X with intensity 

Am{dx) = f3m{x)iK{ym,x) An)n{dx), 
where (3^ : ^ — )■ [0, 1] is given as /3i = 1 and, for m > 1, 



, K(y,x) 
1 _ ^ Al 



n 



, X e X. 



Consider first the case m = 1. Denote the Poisson point process on {X,T) used in the construction of 
RG„(k, /i) by Nn{K,n). Prom the complete independence property of Poisson point processes and the 
non-atomic assumption on /u, conditioned on the existence of a vertex xq in Nn{K,fi), NniK,fi) \ {xq} is 
once again a Poisson point process with intensity n • n{dx) on X. Also, conditioned on Nn{K,n), a given 
type X vertex in Nn{K, fi) would be connected to xq with probability (k(xo, x)/n) A 1. Thus the neighbors 
of Xq, namely Ei, define a Poisson point process with intensity {k{xo,x) An)fi{dx). This proves the above 
statement on with m = 1. 

Consider now m > 1. Since is non-atomic and Um-i^Rm-i consists of only finitely many elements, it 
follows that conditioned on vertices in Rm-i'JUm-i belonging to Nn{K, fi), Nn{K, n)\{Rm-i^Um-i) is once 
again a Poisson point process on X with intensity n-fi{dx). Note that a vertex x £ Nn{K, ^)\(i?m-iUC/m-i) 
is in Em if and only if x is a neighbor of ym and x is not a neighbor of any vertex in Rm-i- So conditioned 
on C/m-i}, the probability that x is in Em equals 

{K{ym,x)/n A 1) • Uy(zii^_^ [1 - (^(y, x)/n) A 1]. 

From this and the fact that the edges in RGn(K, fi) are placed in a mutually independent fashion, it follows 
that the points in Em, conditioned on {Rm-i, Um-i}, describe a Poisson point process with intensity 

nn{dx) • {K{ym,x)/n A 1) • UyeB.^_^[l - (K{y,x)/n) A 1] 
= ^y£R^-i [1 - {i^{y, x)/n) A 1] • {n{ym,x) A n)fi{dx) 

= Knidx). 

Thus conditioned on {Rm-i-,Um-i}-, Em is a Poisson point process with the claimed intensity. 

Next note that one can carry out an analogous breadth first exploration of BP(xo). Denoting the 
corresponding vertex sets once more by {-Rj, f/j, E'j+ijjeNo we see that conditioned on {Rm-i-,Um-i}-, 
Em is a Poisson point process with intensity K{ym, x)n{dx). 

As < j3m{x) < 1 and k A n < k, we can now construct a coupling between BP(xo) and C^'^(xo) by 
first constructing BP(xo) and then by iteratively rejecting each offspring of type x in Em (and all of its 
descendants) with probability 

^ _ /3m{x){K{ym,x) A n) 

l^iUmj x) 

The lemma is now immediate. ■ 



6.5 Proof of Lemma 16.131 

Assume throughout this subsection, without loss of generality, that 

that K is a symmetric kernel. Define, for A; G N, the kernels k^^^ recursively as follows, and for 

all /c > 1 

K(*''+^)(x,y) = / K'^''\x,u)K{u,y)fi{du). 

Recall that {■^^'^''li^i denotes the /c-th generation of BP(xo) and note that it describes a Poisson point 
process with intensity K^^\xQ,y)fj,{dy). This observation allows us to compute exponential moments of 
the form in the lemmas below. 
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Lemma 6.14 Let g : X ^ R_|_ he a hounded measurable map. Fix 5 > Q and letO < e < log(l + 5)/||5'||oo • 
Then 

Ni 

Eexp{e^9(4^'^)} < exp{eil + 5){JCg){xo)}. 

i=l 

Proof: Fix S, e as in the statement of the lemma. By standard formulas for Poisson point processes 

Eexp{e^5(d'^)} = exp{ / k(xo, ^x)(e^^(") - l)iiidu)} 
i=i •'■^ 

< exp{ / k{xo, n)(l + S)eg{u)iJ.{du)} 
Jx 

= exp{e(l + (5)(/Cff)(3;o)}, 

where the middle inequality follows on noting that e^^(") — 1 < (1 + 6)eg{u), whenever eg{u) < log(l + 6). 
m 

Using the above lemma and a recursive argument, we obtain the following result. Recall that Gk = 
Yli^i ^iCi'^^ ) denoted the volume of generation A; where volume is measured using the function 0. 

Lemma 6.15 Fix A; G N and 6 > 0. Given a weight function (p, define (po = (p + + 6yiQ(j). Then 

k k 

Eexp{e^Gi} < exp{e[cP{xo) + ^{1 + SyV^{xo)]} = exp{e<^o(xo)}. (6.19) 

i=0 i=l 

Proof: Define {^jj^^o using a backward recursion, as follows. Let (pk = (p- < i < k 

(pi = (p + {l + d)IC(pi+i. 
Let Ti = ^ = 1) ■ ■ ■ We will show recursively, as I goes from A; to 0, that 

k Ni 



E[exp{6^GJ|J-;] < exp{e^</.Kef)}}- (6-20) 



.(0 

i=l 1=1 

The lemma is then immediate on setting Z = in the above equation 
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When I = k, (|6.20p is in fact an equality, and so (|6.20p holds trivially for k. Suppose now that (j6.20p 
is true for / + 1, for some / E {0, 1, • • • , fc — 1}. Then 

k k 

E[exp{eJ^Ga|Ji] =exp{eGaE[E[exp{e ^ Gi}\Ti+i]\Ti] 



i=l i=l+l 

Ni+i 

< exp{eGaE[exp{6 

i=l 

Ni 

< exp{eGi} exp{e(l + 5) ^ /C0z+i (^f ^)} 

i=l 

Ni Ni 

= exp{6 </)(4« )} exp{6(l + 5)J^ /C0,+i (C?)} 

i=l i=l 
Ni 



1=1 

Ni 

1=1 



For the first inequality above we have used the fact that by assumption (I6.20p holds for / + 1 and for the 
second inequality we have applied Lemma [6.141 along with the observation that e||</>i||oo < log(l + 5) holds 
for all I = 1,2, k, since for all I, cpi < 0o and e||</>o||co < log(l + 6). 
This completes the recursion and the result follows. ■ 

(k) (k) 

To emphasize that cpQ in the above lemma depends on 5 and k, write (po = (p^ . Note that <pg is 
increasing in k. Let (p^ = limfc_j.oo The following corollary follows on sending /c — )• oo in (|6.19p . 

Corollary 6.16 Fix 6 > and e G (O,log(l + (5)/||0^||oo). Then 

E{expeG} < exp{e(l)}{xo)}. (6.21) 

Lemma 6.17 For n G N and x G X 

fC<j){x) < p^'^WfxhHh, where fx{-) = •) and p = p{k). 

Proof: Note that 

/CXX) = / A^(")(x,^x)0(^xHd^x) < II / /,(n)K("-l)(7/,0M^^^)||2||</'||2 

Jx J X 

= ||/C-Vx||2||0||2</5"^'||/x||2||./'||2. 



Now we can finish the proof of Lemma 16.131 
Proof of Lemma EIS Observing that ||(/)||2 < ||(/)||oo^('^')^/^ and ||/^.||2 < ||K||oo/i(A')i/^ we have for 
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6 £ (0, oo) such that {1 + 5)p < 1, and x £ X 

oo 

0Kx) = </>(x) + ^(i+5r/cv(x) 



1=1 



+ ii/.ii2ii<^ii2(X;(i+w~': 



i=l 



where the first inequahty above fohows from Lemma 16.171 Setting 6 = ^, we see 

(1 + 6)p = (1 + A/2)(l - A) < 1 - A/2. 
Using this and that A < 1, we have 

'^K^)<Woo(i + ^^%^)-rfi. 

Let e = log(l + 5)/(2di). Clearly e G (0, log(l + (5)/||</>5||oo)- Using Corollary EH] we now have that 

F{G > m} < exp{— em} exp{e(/>|(xo)} 

< exp{— em} exp{ — ^ — ^} 

< 2exp{-!25|±i)„}. 

Finally, noting that log(l + 5) > |, we have 



log(l + 6) ^ A2 



2di - 8||(/>||oo(1 + 3||k||oo//W)' 
The result follows. ■ 

6.6 Proof of Lemma 16.91 

We begin with a general result for integral operators on general type spaces. 

Lemma 6.18 Let v^p he two mutually absolutely continuous finite measures on a measure space {X,T). 
Let g = dv/dp. Let k : X x X ^ he a kernel. Define another kernel k' : X x X ^ as 



k{x, y) = Y x,y e X. 

Denote by fC [resp. }C' ] the integral operator associated with k [resp. k! ] on L'^ {X ,T , v) [resp. L'^{X,T, p)]. 
Then \\}C\\l2(u) = \\^'\\L2{^l)■ 

Proof: Note that the operator A : L'^{X,T, v) L'^{X,T,p) defined as {Af) = ^f, f £ L'^{v), is 
an isometry. Also, for / G L'^{fi) 

{A1CA~^ f){x) = y/g{x) i K{x,y)—l=f{y)p{dy)= k {x,y)f{y)fi{dy). 
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Thus AICA"^ = K'. The resuh now fohows on noting that ||^/C^ II-L^(m) ~ ■ 

For the rest of this subsection we wih take {X, T) to be the cluster space introduced in Section 15.41 
(see above ()5.17p ). Given rate functions {a,b,c), fit{a,b,c), Kt{a,b,c), (j)t{a,b,c),ICt{a,b,c), pt{a,b,c) areas 
introduced above Lemma [6191 



Lemma 6.19 Let {ai,bi,Ci), i = 1,2 be two sets of rate functions. Suppose thatai,ci are strictly positive 
on (0, r] and 02 < ai,C2 < ci on [0,T]. Also suppose that for some 6 G (0,e~^/r), ci < C2 + S on [0,T]. 
Fix t G [0,T]. Let fj.i = fj.t{ai,bi,Ci), i = 1,2. Then <C fJ-i and 



^{s,w) = ^ X exp{- [ w{u)[c2{u)-c,{u)]du} x ngP~' (^) , {s,w) G [0,t] x W. 
where Ti = Ti{s,w) is (fii a.s.) the i*^ jump of w after time s. 

Proof: Recah the probabihty measure on W introduced below (|5.17p . Write vl = Usia-i, bi,Ci), i = 1,2. 
Note that (see (|5.18p ). for s G [0, t], m^ds dw) = vl{dw)ai{s)ds, i = 1,2. Thus to prove the result it suffices 
to show that, for all s G [0, t], i>1 <^ and ^ = LJ, where, for t G [s, T], 

L\{w) = exp{-^ w{u)[c2{u) - ci{u)\du] x nj>i ^^^|^l{a,<t}^ , 

and (Ti{w) is {vl a.s.) the i*^ jump of w. For this, it suffices in turn to show that jy^ L\{w)vl{dw) = 1 
for all t G [s,T] (see eg. Theorem T3, p. 166 of [8j). The process {L\}t£[s,T\ 011 ^ with the canonical 
filtration is a local martingale under (see Theorem T2, p. 166, [8J) so to check the martingale property, 
it suffices to check (see (2.4) on page 166, and Theorem T8 on page 27 of [8]), that 



I [[ L'^s{w)\ci{u) - C2{u)\du]vl{dw) < +00. 

Jw Js 



Note that L^{w) < ex.p{Tw{T)5} since ci — 5 < C2 < ci. Thus 



L'^{w)\ci{u) - C2{u)\du]ul{dw) < 2T / e^'"^^'^^u^{dw). 



Note that, under Uq, w{T) is stochastically bounded by the sum of two independent Geom(e~^) (see 
proof of Lemma l6.7p . Thus the last integral is bounded by (Ee'^^^)'^, where Z is a Geom(e~"^) random 
variable. This expectation is finite since T5 < e~^. The result follows. ■ 

Proof of Lemma [6.91 will make use of parts (iv) and (v) of Lemma [6.51 In order to use (v) we will need 
the kernels to be suitably bounded. For that we use a truncation of the kernels, the associated error of 
which is estimated through the following lemma. 

Lemma 6.20 For A G (0, 00), t G [0, T] and rate functions (a, b, c), define the kernel KA,t(,CL, b, c) = KA,t '■ 
X X X ^ [0, 00) as 

KA,i(x,y) = Kt(x,y)l{^(T)<A,tD(T)<A}, where Kt = Kt{a,b,c),yi = {s,w),y = {r,w). 

Denote by K,A,t the integral operator corresponding to K,A,t on L'^{X ,T, fJ-t) CLnd p{KA,t) its norm, where 
pt = Ptia,b,c). Then there exist Ao,C3,C4, G (0,oo) such that 

p{Kt) - Cse"^*^ < p{KA,t) < pint), 

for all rate functions {a,b,c), t G [0,T] and A > Aq. 
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Proof: We will suppress t in the notation. Since ka < k, from Lemma 16.51 (ii) and (iii) we have 

p{k) - p{k - ka) < p{ka) < 

Consequently, 

p{k) - p{ka) < p{k - ka) < ||k - KAh, 
where ||k||2 denotes the L^(/U x p) norm of the kernel k. Note that for x,y £ X 

«;(x, y)= / w{u)w{u)b{u)du < Tw{T)w(T), for p x p a.e. (x, y) = ((s, z/;), (r, tD)). (6.22) 
Jo 

Let A = {{s, w) e X : w{T) > A}. Then for fixed x = [s, w) e 

{K{x,y) - KAix,y)fp{dy) < / {K{x,y)lA{y)f p{dy) 

X J X 



< / [TW{T)w\T)l{^i^T)>AMdrdw)) 
Jx 

< tW{T)Eo[w^{T)1^^^t)>a}], 



where in the second line we have used (j6.22|) and Eq in the third line denotes the expectation corresponding 
to the probability measure z^o on W. Next, for fixed x = (s, w) £ A, we have in a similar manner 

/ (K(x,y) - KA{x,y)fp{dy) < T^w\T)Eo[w\T)]. 
Jx 

Thus for any x = (s, w) G , 

/ (K(x,y) - KAi^,y)fp{dy) < TW{T)Eo[w\T)1{^^t)>A}] + l{y,iT)>A}T^w^{T)Eo[w\T)]. 
Jx 

Integrating with respect to x G Af, and noting that i's{w{T) > a) < i'q{w{T) > a), a > 0, we have 

/ / ii^i^,y) - tiA{y^,y)fpidy)pidyL) 
Jx Jx 

<T''Eo[w\T)l{^^T)>A}] [ a{s)Eo[w\T)]ds + T^Eo[w''iT)] [ a{s)Eo[l{^^T)>A}w\T)]ds 

Jo Jo 



{k - ka)\\1 



< 



2T'Eo[w\T)]Eo[w\T)l^^^T)>A}]. 



As noted in the proof of Lemma [6.71 under fg, w{T) is stochastically dominated by + Z| where Z*, Z| 
are two independent copies of Geom(e~-^). Therefore 

Eo[wl{T)] < E{Zl + Z^f = 26^^(2 - e"^) < 46^^ 

and for a suitable Aq £ (0, oo) 

IEoK(T)l{.„(T)>A}] <Y.ex{k- l)e-2^(l - e-^r^ < (1 - 6"^)^/^ 

k>A 

for all A > Aq. Combining the above estimates, for all A > Aq 

||Ac-ACA||2<8rV^(l-e-^)^/2^ 

The result follows. ■ 

In the proof of Lemma 16.91 we will apply Lemma 16.181 to measures pi,p2 of the form in Lemma 16.191 
This latter lemma shows that (under the conditions of the lemma) p2 <^ pi- However, to use Lemma 
16.181 we need the two measures to be mutually absolutely continuous. To treat this difficulty we will use 
an additional truncation introduced in the lemma below and the elementary fact in Lemma 16.221 
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Lemma 6.21 Given rate functions {a,b,c) and A G (0,oo) and t G [0, T], let Kt^KA,t be as in Lemma 
\6.2(h For {s,w) £ X, let t{s,w) = inf{n > s : w{u) - w{u-) / 0}. For 5 > 0, let As = {(s^w) £ X : s < 
S and t{s,w) < s + 5}. Define the kernel KA,s,t cls 

KA,5,t (x, y ) = KA,t (x, y ) 1ai (x) 1ai (y ) • 
Then there exists a £ (0, oo) such that for all rate functions (a, b, c) and 5,A£ (0, oo) 

PiKA,t) - C^A^S < p{KA,5,t) < PiKA,t)- 
Proof: Once more we will suppress t from the notation. As in the proof of Lemma 16.2^ we have 

p{ka) - \\ka - KA^sh < p{ka,s) < p{ha)- 

For X G A"", 

\2 . ^ I .,2 

I- 
10 



(KA(x,y) - KA,5(x,y)) ^(dy) < / KA{yi,{r,w))l{r<r+5}'^{T<&}P{drdw) 
X Jx 



< T^A"^ / Ur{T{r, w) <r + 6}dr, 
Jo 



where the above inequality uses the bound ka < TA^. 
Also, for fixed x G A, 

/ (KA(x,y) - KA,5(x,y))V(rfy) < T^A^ ■ T 
Jx 

Thus for all X G A', 

/ (ka(x, y) - KA,six, y)) V(t^y) < T^A* I ^r{T{r, w)<r + 5]dr + lA(x)r3^^ 
Jx Jo 

Finally 

\\'^A-KA,s\\l= / (KA(x,y) - KA,5(x,y))^^((iy)^((ix) < 2r^A^ / Ur{r{r,w)<r + S}dr. 
Jx Jx Jo 

The result follows on observing that i'r{T(r, w) < r + 6} = 1 — exp{— 2c{u)du} < 26. m 

We will use the following elementary lemma. 

Lemma 6.22 Let 70, 71 be finite measures on a measure space {X,T) such that 70 <^ 71. Let G £ T be 
such that {x : djo/dji > 0} D G. For i = 0,1, let 7^ be restriction ofji to G: 7p(-) = ji{-nG). Then j^f 

and 7p are mutually absolutely continuous with ^p-(x) = ^(x)Ig'(x) and ^^(x) = (^(x))^"'^1g'(x), 

a.s. 7^* and . 

Lemma 16.91 requires establishing an estimate of the form C2(— log (5)^(5^/^ for both \pt — p+,t\ and 
\pt — P-,t\- Proofs for the two cases are similar and so we only provide details for \pt — P-^t\ and leave the 
other case for the reader. Given 5 £ (0, 00) and rate functions (a, b, c) as in the statement of Lemma 16.91 
we denote {as, bs, cs) = ((a — d)~^, {b — 5)+, (c — 5)"*"). Note that 

\pt-P-,t\ = \pt{a,b,c) - pt{as,bs,cs)\ 

< \pt{a,b,c) - pt{as,b,cs)\ + \pt{as,b,cs) - ptias,bs,cs)\. (6.23) 

We treat the first term on the right side in Lemma 16.231 while the second term is estimated in Lemma 

MM 
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Lemma 6.23 Let (a, b, c) be rate functions as in the statement of Lemma \6.9l There exists Cg, 62 G (0, 00) 
such that for all 6 G (0, 62) and t £ [0, T] 



\Pt 



(a, b, c) - ptias, b, cs)\ < Cei- log 6f5^/\ 



Proof: Since the kernel Kt{a,b,c) does not depend on a,c, nt{a^b,c) = Kt{as,b,cs) = Kt- Henceforth 
suppress t from the notation. Let p = pt{a, b, c) and ps = pt{as, b, cs). For e > 0, let C be defined as 

D = De = {{s,w) : w{T) > A or (r(s, w) < s + e and s < e)}" 

and define the kernel as 

K'D{^,y) = K(x,y)lz)(x)lD(y). 
Using Lemmas 16.201 and Lemma I6.2H 

p - Cse-^^^ - C^Ah < p{fiD) < P 

and 



PS 



g{s,w) =a,f ^{s,w) = ^^44 X exp{- / w{u)[cs{u) - c{u)]du} x II^^P ^ ( - 
dp ao[s) Jg \( 



where p{hd) [resp. ^^(k/))] is the norm of the corresponding integral operator on L?'{p) [resp. L?'{ps)\, 
where p = p{a, 6, c) and ps = p{as, 6, cs). Thus we have 

\p - ps\ < 2C3e-^^^ + 2C5Ah + \p{kd) - Ps{kd)\. (6.24) 

We now estimate \p{kd) — psi'^D)\- By Lemma [6. 191 ps ^ P and 

p a.s., where are as in the statement of Lemma 16.191 For p a.e. {s,w) £ D we have 

w{t) < w{T) < A and ri(s, w) > e, 

consequently cs^Ti) > {9e — 5)~^ for all i, where 6 is as in the statement of Lemma 16.91 Also, since a is 
bounded away from 0, we can find di G (0,oo) such that as{s) > {di — 5)'^. Thus for p a.e. x £ D 

{di - 5)+ {{Oe - (5)+)^ < 5(x) < exp{TA6}. (6.25) 

Denote by p^ [resp. p^] the restrictions of p [resp. ps] to D. Then from Lemma 16.221 whenever 
S < 60 = mm{9e,di}, p^ and pf are mutually absolutely continuous and 

dp^ 

-^(x) = fif(x)lz)(x), a.e. dp^. 

For the rest of the proof we consider only 6 < 5q. Then, 

5 . 5^^ 



(1--)(^1--J <gix) <exp{TA6}. (6.26) 

Note that p{K£),p) = p{kd,P^) and p{kd,Ps) = piKDjpf). Also by Lemma [6.18|, p{K£),p^) = 
p{k'j^, p^), where 

«^D(x,y) = KD{^,y)J ^|^lD(x)lz)(y), x,y G Af. 
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Thus using Lemma 16.51 



\p{f^D) - Ps{i^d)\ = \p{KD,fJ.) - piKD,P5)\ 

= \p{'^D,P^) - p{kd,Ps)\ 

= \p{l^D,f^^) - p{K'r),P^)\ 



9[y) 



< T sup ( ij - l|Kz)(x,y) 



9iy) 



Prom (|6.26p we see that whenever 6 < 



sup ( V I < (1 + ^6) X eMTA6} x (1 + f)^ ^ d{6,e,A). 



Noting that kd < TA^, we have 



sup - l|KD(x,y) < TA^\d{6,e,A) - 1|. 

Thus 

Combining this with (j6.24p . we have 

\p-ps\ <2C3e-^^^ + 2C5A^e + TA'^\d{6,e,A) (6.27) 

Take A = — ^log(5 and e = J^/^. Note that when 5 is sufficiently smah, 5 < ^minj^e, di} and so the 
above inequahty holds for such 6. Also, with this choice, we can find 6i,d2 G (0, oo) such that the sum of 
the first two terms on the right side in ()6.27p is bounded by 

2C3(5 + ^(- log (5)25^/2 < d2i-log6f6^/^ for all 5 < di. 

Also note that with the above choice of £,A, d{5,e,A) — )■ 1 as (5 — t- 0. Furthermore 

di5,e,A) = (1 + 0{5))il + 0{5)){1 + 0{5^/^i- log<5))). 

Thus,we can find d3,62 G (0,oo) such that whenever 6 < 62 

TA'^\d{6,e,A) - 1| < d3{-log6f6^/^. 

The result follows on combining the above estimates. ■ 
We now estimate the second term in (j6.23p . 

Lemma 6.24 There exists C7 G (0, 00) and 63 G (0, 1) such that for all t G [0, T] 

\pt{as,b,cs) - ptias,bs,cs)\ < C7(- log 5)^5 



34 



Proof: Denote Kt{as,bs,cs) = ns^t = ^5 and recall that Kt{as,bs,cs) = Kt = k. We suppress t in rest of 
the proof. For A > 0, let ka [resp. ks^a] be the truncated kernels as in Lemma 16.201 associated with k 
[resp. Kg]. Clearly 

\l^A — 1^5, a\ < TA'^5. 
Using this bound and Lemma 16.201 for all A > Aq 

\p{as,b,cs) - p{as,bs,cs)\ = \p{k, p^s) - pi^s, ns)\ 

< TA^5 + \p{k, ps) - p{i<'A,P5)\ + \p{i^&, Ps) - p{ks,a, Ps)\ 

< TA^6 + 2Cse^p{-C4A}, 

where ps = pi^s, bs,cs) = p{as, b, cs) is as in Lemma [6.231 

Result follows on taking A = (— log(5)/C4 and S sufficiently small. ■ 

Proof of Lemma 16.91 The proof is immediate on using Lemmas 16.231 and 16.241 in (j6.23p . ■ 
6.7 Proof of Lemma 16.101 

We begin with an elementary lemma which allows one to regard the operators /Cj (a, b,c), t G [0, T], to be 
defined on a common Hilbert space. Recall that for a kernel k on X x X and a finite measure p on X , 
we denote by p{k, p) the norm of the integral operator associated with k on L'^(X ,T, p) = Lp'{p). 

Lemma 6.25 Let (a, 6, c) be rate functions. Then for all t € [0, T] 

p{Kt{a,b,c),pt{a,b,c)) = p{Kt{a,b,c), pT{a,b,c)). 

Proof: Write Kt{a,b,c) = Kt, Pt{a',b,c) = pt- Denote the integral operator corresponding to on L^{pt) 
[resp. L'^{pt)] by ICt [resp. ICf]. Let Xt = [0,t] x W, then pt is supported on Xt and Kt is supported on 
Xt X Xt. Thus for any -0 G L,'^{pt), f^fi^ = is also supported on Xt, and this implies \\}Cf\\ < \\JCt\\. 
On the other hand, for any -0 £ L'^iPt), = ^t"(^l[o,i]xw) 

\\ICt^\\mi,t) < ll^fllllV'l[0,t]xwllL2(/.T) = ll^fllllV'llL2(/.t)- 

Thus ||/Cf II > \\ICt\\. m 

The following theorem concerning IRG models on a general type space is a corollary of Theorem 3.1 
and Theorem 3.12 in |16| . 

Theorem 6.26 116^ Let {X,T,p) be a type space. Consider the weight function (j) = 1 on this space. Let 
Kn, K be symmetric kernels on X x X such that 

X and Vn ^ U implies K„(x„,y„) K{x,y). 
Let Cn\K) denote the size of the largest component in'RGn{K, p,<j)). Then 

(i) If p{k, p) < 1, then Cn {ku) / n — >0. Furthermore, i/ ||k||oo < +00, then Cn {ku) = O {log n). 

(ii) If p{k, p) > I, then Cn\Kn) = 6(n). 

Proof of Lemma 16. lot Let Kj(ao, ^Oi cq) = Kt and pt{cLo, bQ,co) = pt. Note that for < ti < ^2 ^ 
since Kt^ < Kt^, we have from Lemma 16.251 that 

P(.h) = p{Kti,Pti) = p{Kt^,PT) < p{Kt2,PT) = p{Kt2,Pt2) = P{t2)- 
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Thus p{t) is nondecreasing in t. Next note that, since w{u)w{u)bo{u) is non decreasing in u for fiT x fJ-T 
a.e. ((s, w), (r, tD)), Kt(x,y) is convex in t for a.e. x,y, i.e. for ht x a.e. (x,y) G X x X, and all 

ti,t2 G [0,r], and a,/3 G [0,1], Q + /3 = 1, 

«^ati+/3t2(x,y) < aKf,(x,y) +/3Kt2(x,y). 

Thus 

p(Qtl + /3t2) = p{Kati+l3t2,fJ'T) 

< p{atiti +/3Kt2,/iT) 

< p{aKt-^ ,Pt) + p{fint2 , /^t) 
= ap(Kt, , //t) + Pp{i^t2 , /^t) 
= ap{ti) + /3p(t2), 

where lines 3 and 4 above use parts (ii) and (iii) of Lemma 16.51 and line 2 uses the convexity of k.. 

Thus p is convex on [0,T]. Also since p(0) = and p{t) > for t > 0, we have that p is strictly 
increasing on [0, T] and has a strictly positive left derivative on (0,T]. This proves parts (i) and (iii) of 
Lemma 16.101 

We now consider part (ii). For 6 > 0, let 

= Ptiiao + 5) A 1, ibo + 6)M, (Co + 5) A 1), /-(t) = ptHao - 5)+, (6o - <5)+, (cq - <5)+). 
Similarly define fif~, i^f~, P'f'^ and k^'^. We will argue by contradiction. 

Suppose first that pitc) > 1, then by Lemma [6.91 and the continuiity of pit), there exist e,5 > 
such that p^~itc — e) > 1. Denote RG^'~(k„) = RG„(k^'7, </>*), where is defined as in (I5.19p . 
replacing b there with (6o — <J)^- Since converges to uniformly on compact subsets oi X x. X , by 
Theorem 16. 26|, the size of the largest component of RG^' (^n); whp, is 0(n) and consequently the volume 
of the largest component of RGj'~(k„) is, whp, at least G(n). By Lemma and as x„(t) — >• x(t), we 
have whp 

COMnitc - e) >d IA„((ao - 5)+, (6o - 5)+, (co - 5)+)t.-. 

and from Lemma [5.31 the largest component in IA„((ao — (5)"'",(6o — ^)''')(co — ^)~^)tc-e has the same 
distribution as that in 'RG^^'^^inn). However, by Theorem 1.1 of |41] the largest component size of 
Bohman-Frieze model for t < tc is, whp, 0(logn) which contradicts the fact that the volume of the 
largest component of RG^'7e(K„) is, whp, at least 0(n). Thus we have shown that pitc) < 1. 

Suppose now that pitc) < 1- then there exists e,6 > such that p^^'^itc + e) < 1. Then a similar 
argument as above shows that, whp, 

COMnitc + e) <d IA„((ao + <5) A 1, (5o + <5) A 1, (cq + 5) A l))t^+, =def IA^,t,+e- 

Lemma [6.111 implies that whp the largest component in IA^^^^_^^ is O(log^n). However from Theorem 1.1 
of [S], for t > tc, the largest component size of Bohman-Frieze model is whp 0(n). This contradiction 
shows that pitc) > 1. Combining the above arguments we have pitc) = 1. ■ 

7 Proof of Proposition 14.1! 

We shall now study the asymptotics of ^2 and ^3, namely the sum of squares and cubes of the component 
sizes. We first analyze ^2 since the asymptotics for this will be required in the analysis of ^3. 
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7.1 Analysis of S2( ) near criticality 

Let us start with the sum of squares. Recall from (|2.1|) that 52 (i) denotes the sum of squares of the 
component sizes in BF(t) and S2{t) = S2{t)/n. Also recall the limiting functions Sk{t), k = 2,3, introduced 
in (j2.3p and (j2.4p . Note that these functions are non-decreasing and they blow up at the critical point tc 
(see (j2.7p and ()2.8p ). Let y{t) = l/s2{t). Then this function satisfies the differential equation 



y'{t) = -x\t)y\t) - (1 - x\t)), y{0) = 1, t G [0, t^]. 



(7.1) 



Note that y is a monotonically decreasing function with y{t) — t- as i — t- and as shown in Theorem 3.2 
of [26], the scaling behavior near tc is 



1 



y{t) = -{tc-t) + 0{tc-t) 



a 



(7.2) 



as t t tc, where a is as in ()2.5p . Let Yn{t) = l/s2{t). To simplify notation, we suppress the dependence 
of the process y on n when convenient. Note that ()4.2p is equivalent to showing 



n 



1/3 



0. 



Here, and throughout Section[71 tn = tc — l/n"^ with 7 G (1/6, 1/5). From ()7.2 



y{tn) 



1 



Thus to show (j7.3p it is enough to prove the following 
Proposition 7.1 As n ^ 00 



n 



sup \Y{s)-y{s)\ 



n 



0. 



1/3 



(7.3) 



(7.4) 



s<t„ 



We shall prove this via a sequence of lemmas. We begin with some notation. Recall that C^^(t) = 
{Cn\t) : i > 1) = {Ci{t) : i > 1) is the component size vector, In{t) the size of the largest component, 
and Xn{t) the number of singletons, in BF„(t). Let Yli denote the summation over all components and 
denote the summation over all pairs of components with i < j. 

The first Lemma identifies the semimartingale decomposition of the process Y(-) as well as the pre- 
dictable quadratic variation (M) of the martingale M in the decomposition. Recall the natural filtration 
associated with the BF process introduced in Section | 



Lemma 7.2 The process Yn{-) can be decomposed as 



Yn{t) = l+ [ An{s)ds + Mnit), te[0,tc 

Jo 



where 

(a) Mn is a RCLL martingale with respect to the natural filtration {J-t\t>Q of the BF process, 
(h ) The process An = A^ + i?" where ( suppressing n ) 



Ai{u) = -y2(n)x2(n) - (1 - x2(n)) + (1 - x2(u)) 



y2(n) 



^.4 

i 



Cf{u), U<tc 



and for some C7 G (0, 00), 



\RUu)\ < Cj 



1 I IrXu) 

n n 



for all n and u < tc- 



(7.5) 



(7.6) 



(7.7) 
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(c) Predictable quadratic variation of Mn is given as 



where Bn is such that 



{Mn){t) = / Bn{u)du 

Jo 



Bn{u) < - + ' ^ U < t,. 

n n 



Proof: We will suppress n from the notation when convenient. Note that 

Y{t) = l + Y^ Ay(s), where AY{s) = Y{s) - Y{s-). 



s<t 



We now analyze the possible jumps of Y . Note than any jump in Y corresponds to a jump of one of the 
Poisson processes "Pe, e = (61,62) G 1?^ (recall the notation from Section [2D. A jump of Ve at a time 
instant u could result in two different kinds of jumps in Y . 

(i) Merger caused by the first edge ei: In this case /S.S2{u) = S2{u) — S2{u—) = 2 which implies 



Ay(M) = ai{u- 



2Y'^{u-) 



n 



1-0 



2Y{u-) 



n 



(ii) Merger caused by the second edge 62: In this case, suppose components i and j merge, then 
AS2{u) = 2Ci{u—)Cj{u—) and thus 



n 



1-0 2 



Ci{u-)Cj{u-)Yniu-) 



n 



With these observations we can represent Y in terms of stochastic integrals with respect to Ve as follows. 
Define 

7ii{u) = {e = (61,62) G f"^ : 61 = (^1,^2) where both 'i;i,f2 are singletons at time u}, 



nf^\u) = {e = (61,62) G £\Hi{u) '■ 62 = {vi.,V2) where one vertex is in C^{u) while the other is in C-'{u)^. 
Also let 



(n) = ai(^/)l^^(„)(e), Ul;^ {u) = aWu)l^(,,,)^^^^{ 



e . 



Then 



y{t) = l+Y. I (^e(.-)+E^e^'(^-))^e(rf^). 



Recalling that Ve is a rate Poisson process, one can write Y as 



where 



Y{t) = 1 + / A{s)ds + M{t), 



(7.9) 
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Note that 



-xl{s)-{l + 0{l/n)), 



(7.10) 



2 



n 



C^(s)C,(s) = y (1 - xlis))Ciis)C,{s) ■ (1 + 0(l/n)).(7.11) 



Thus we get ^(s) = Ai{s) + Ri{s), where Ai represents the leading order terms: 



n 



2y2(s) 



Ms) = -i,xi{s) ■ 

2 n 



1 



Y^-{l-xl{s))Ci{s)C,{s) 



2Y^{s)Ci{s)Cj{s) 



i<j 



n 



-X. 



-x: 



,)Y\s) - il-xlis))Y\s) . -L^2Cf(s)Cj(s) 

i<j 

l{s)Y\s) - (1 - xl{s))Y\s) • ^[{Y.C!{s)f - Y^Cfis)] 

i i 

s)Y\s) - (1 - xlis)) + (1 - xl{s))Y\s) • EC'(^) 



and the last line follows from the fact Y 



The term Ri consists of the lower order terms and 



using the observations Y < 1, x^ < I and < x'^Y'^ + (1 ~ x'^) ^ 2, it can be estimated as follows. 



in / M ^ I ^ / M di n_n, , 2y^(n) d2Y(u) 

\Ri{u)\ < \Ai{u)\ ■ + -xliu) ^ • -^-^ 

n 2 n n 

+ (1 - xl{u))Y\u) ■ \ Y.'(^CKu)C]{u) ■ 



dj,Ci{u)Cj{u)Y 



■' — ' n 

Kj 



u)] 



n n 



n 



i<j 



n n 

Next, from (|7.9p and using independence of Poisson processes Ve, 



B{s)ds, 



where, using (j7.10p and (j7.1ip once more, we can estimate B as follows. 



no,.. l2Y^{u)\^ ^ ^2,.,^^ o ^Q{u)Cj{u) [2Y^{u)Ci{u)Cj{uy'^ 



B{u)<-xi{u)-2 

2 \ n 



+ {l-xl{u))-2.Y, 



4 AY^{u)P{u) Y^iu) ^ N^2/ N 

Kj 



i<j 



n 



n 



^ 4 4Y^{u)P{u) 
~ n n 

This completes the proof of the lemma. 
■ 

The following result bounds the difference Y — y through an application of Gronwall's lemma. 
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Lemma 7.3 For every n € N, 



sup \Yn{s) - y{s)\ < e„e 

s<t„ 



2tc 



where 



En = 4:tc sup \x{s) - x{s)\ + SUp |M„(t)| + Cg 

s<tc s<tn Jo n 



Cg = (2C7 + 1) and Cj, Mn are as in Lemma 1.2 



Proof: Since y solves ()7.ip and Y satisfies (j7.5p . using the fact tliat x,x take values in [0, 1], we have for 
any fixed t < tn 



\Y{t)-y{t)\< 



{Y\s)x^is) - y^is)x\s)) ds 



t r-t 

+ / \x'^{s)-x'^{s)\ds+ / \Ri{s)\ds 
Jo 



/O \ "-" ~ / s<t 

Let us now analyze each term individually. Writing 

Y\s)x\s) - y\s)x\s) = xHs){Y\s) - y\s)) + yHs){x\s) - x^s)) 
and using the fact that Y, y, x, x take values in the interval [0, 1] we get 



Y^{s)x'^{s) -y^is)x'^{s)ds 



<2 [ \Y{s) - y{s)\ds + 2tc sup \x{s) - x{s) 

Jo S<tc 



Similarly the second term in ()7.12p can be estimated as 

/ \x'^{s) — x'^{s)\ds < 2tc sup \x{s) — x{s) 

Jo S<tc 

Prom (|7.7p the integrand in the third term in (j7.12p can be bounded as 

|;j,(,)|<C,fi+«£ll<2c/"« 



n n J n 
Similarly, the integrand in the fourth term in (|7.12p can be bounded as 

'''^'^Y.Ct{s)<'MY\s)lY.CKs). 



n 



Noting that by definition Y[s)\ YliCfis) = 1 and Y{s) < 1 we get 



n 



Combining we get, for all t < tn 



\Y{t)-y{t)\<en+ f2\Y{s)-y{s)\ds. 
Jo 



(7.12) 



Thus Gronwall's lemma (see Theorem 5.1, Appendix in [23]) proves the result. ■ 

One last ingredient in the proof of Proposition 17.11 is the following estimate on the martingale in 
Lemma 17. 2i 
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Lemma 7.4 As n ^ oo, 

sup \M{s)\ 0. 

s<t„ 

Proof of Lemma 17.41 is given in Section 17.21 



Proof of Proposition I7.lt In view of Lemma 17.31 it is enough to show n^^^En as n — t- oo, 
where En is as defined in the Lemma l7.3[ Let us analyze each of the terms in En- First note that by 
Lemma 16.41 for any i!) < 1/2, and in particular for -i? = 1/3 

sup n'^|x(s) — x(s)| 0, as n — )• oo. 

S<tc 

Next, from Proposition 15.11 

{/„,(t) < m{n,t),\/t < tc — n~'^} whp as n — )• oo, 
where m{n,t) = B ■ Thus, recalling that 7 G (1/6, 1/5), we have, whp. 



on n Jo {tc- s) 

n \n^l-^ J 

The result now follows on combining these estimates with Lemma 17.41 ■ 

7.2 Proof of Lemma 17.41 

We shall prove the lemma by first showing that n''sup^<^^ |M(s)| -^-^ for any ■!? < 1/5 and then 
sharpening the estimates to allow for any < 1/3 . 

Lemma 7.5 As n ^ 00, 

^,,p|M..(.)|>=!Z«Uo (7.14) 



and 

^Ynit) < 2y{t) yt<tn)^l. 
Proof: Fix 71 G (1/3, 1/2) and define stopping times Tj, i = 1, 2 by 

n = inf{t : I^(t) > m{n,t)}, T2 = inf{t : \x{t) - x(t)\ > n^'^^}. 

From Proposition 15.11 and Lemma [67 



P(ti A T2 > t„) ^ 1 (7.15) 

as n — )• 00. Let r* = tn A ri A T2. For (j7.14p . in view of (|7.15p . it suffices to prove the statement with t„ 
replaced by r*. By Doob's maximal inequality we have 



Furthermore, from (17.8 



E(sup |M(s)p) < 4E((M)(t*)) = 4E / B{u)du. 

S<T* Jo 



n n n n 
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Since for t < r* we have In{t) < B{logn)'^ /{tc — t)^ (see (|5.ip ). we have 



Combining the estimates, we have 



E(sup |M(.)|)<dl/^"''^'^^°^")' 



n 



s<T* 

A simple appUcation of Markov's inequahty and (I7.15|) gives (I7.14p . To get the final assertion in the 
Lemma, note that on the set 

f n'^'^'/^flosn)^ 1 

Bn = {n > tn} n {t2 > tn} n <^ sup \M{s)\ < ^ B ; I 



s<T* Vn 



\x{t) - x{t)\ < n-^i and (see fTTSD ) 



I 



Therefore, since 7 < 1/5, the error e„ in Lemma 17.31 satisfies whp 

^37/2 Q N? / I 
£n < V- =0 

/n V nT 



Noting that y{t) is a monotonically decreasing function, the above along with ()7.4p implies that Y{t) 
(1 + Op(l))y(t) for t < i„ on the set Bn- Since P(-Bn) — ^ Ij the result follows. ■ 



Proof of Lemma 17. 4t Along with the stopping times ti,T2 introduced in Lemma 17.51 consider the 
stopping time 

T3 = mf{t:Y{t)>2y{t)} 

Then Lemma 17.51 and ()7.15p imply that 

P(ri A rs A rg > t„) ^ 1 

as n — )• 00. Thus to complete the proof, it is enough to show that for the stopping time r = t„ Ari Ar2 AT3, 
n-^/^E(supg<^ M(t)) — )• as n — )• 00. Once again by Doob's maximal inequality it is enough to show that 

n2/^E((M)(r)) ^ as n ^ 00. (7.16) 

Now note that ()7.8p implies that 



E M r <e/ - + ^ ^ M 

Jo \n n J 



< l^Ii + 1 ^^2(.) 2i32(logn)8 ^_^ 

<4 a+irt"°s">;d.+"°^"> 



In the second line of the above display we have used the fact that /n(t) < m(n,t) for all t < ti and in the 
last line we have used (17. 2p . Thus 

Since 7 < 1/5 we have (|7.16p and this completes the proof. ■ 
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7.3 Analysis of S3 near criticality 

We will now analyze the sum of cubes of component sizes near criticality and consequently prove (j4.1|) . 
Define 

z{t) = ^,t€[0,t,). 

Then differential equations (12. 3j) (12. 4p imply (see [26]) that z solves the differential equation 

z'{t) = 3x\t)y^{t) - 3x\t)y{t)z{t), z{0) = 1, t € [0,tc) 
and furthermore z{t) — )• /3 as t — )• tc- Now consider the process 



Zn{t) 



Sz{t)/n 
iS2{t)/ny 



n 



Then to show (j4.ip . it is enough to show the following proposition: 
Proposition 7.6 Fix any 7 € (1/6, 1/5) and let as before tn = tc — . Then 

\Z{tn)-z{tn)\^^ 



(7.17) 



as n —)■ 00. 



The analysis is similar to that for S2{-) as carried out in the Section 17.11 We begin by writing the 
semimartingale decomposition for Z„ and identifying the predictable quadratic variation (M) of the 
martingale M in the decomposition. 



Lemma 7.7 The process can he decomposed as 



Zn{t) = 1 + / An{s)ds + M„(t), t £ [0,tc] 

Jo 

where 

(a) Mn is a RCLL martingale with respect to the natural filtration {/"j}t>o of the BF process, 
(h) The process — ^4" + i?™ + R2 where (suppressing n), for u G [0, tc], 



ii(n) = 3x'^{u)Y^{u) - 3Z{u)Y{u)x'^iu), 



\Ri(u) 



3{l-x^{u)) 



Z{u)Y{u) 



Cfiu) 



< 



6/3(n)y2(^i) 



n 



and 



~ „ ^ 3QY\u)ll{u) , l2Y\u)I^^{u) 

\R2{u)\ < \ ^ 



n 



(7.18) 



(7.19) 



(c) Predictable quadratic variation of Mn is given as 

{Mn){t)= [ Bn{u)du, 

Jo 

and the process Bn can be bounded as 

B(u) < C, (IIMIM. + yHu)lHu) ^ Y'{u)lHu) ^ Y\u)I^\u) ^ Y\u)I^^{u) \ 

for some Cg G (0,oo). 
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Proof: We will proceed as in the proof of Lemma |7.2[ Note that 

Z{t) = 1 + ^ AZ(s), where AZ{s) = Z{s) - Z{s-). 

s<t 

Next note that if A52(u) = a and lS.S^{u) = b, then 

^ ssju-) + b/n ^ Y^{u-){s3{u-) + b/n) 
(s2(n-) + a/n)3 {I + aY{u-)/n)^ ' 

Using the estimate |(1 + x)~^ — (1 — 3x)| < 6x^, for < x < 1, we get 

Ziu) = Y^u-) (^ss{n-) + ^) (l - ^^^^ + ^(«, , 

where 

\R{a,u-)\ < 

Thus 

A ^/ N , ^ 3aZ(u-)Y{u-) Y^{u-)b ~ . . 

AZ{u) = C(a, b, U-) = ^ ^ + — ^ ^ + Ra,b{u-), 

n n 

where the remainder term 

IP f \\^7( Mfi( ^, , ^abY\u-) , bY\u-)\R{a,u-)\ 
\Ra,b{u-)\ < Z{u-)\R[a,u-)\ H 5 \ 



Qa^Z{u-)Y'^{u-) 3a6y^(n-) 6a26y5( 



< ^ k ^ - + ^ - + 



u- 



Any jmnp in ^2,53 or Z corresponds to a jump of one of the Poisson processes V^, e = (ei, 62) G <?^. A 
jump of Ve at a time instant u could result in the following different values for o, b. 

(i) Merger caused by the first edge ei: In this case a = 2 and 6 = 6 and Ra^b can be estimated as 

,■- , 2AY^(u-)Iiu-) 36y^(u-) UW^{u-) 2MY^iu-)I{u-) 
\R2fi[u-)\ < 2 + 2 + 3 - 2 • 

(ii) Merger caused by the second edge 62: In this case, suppose components i and j merge, then 

a = 0ij{u—) = 2Ci{u—)Cj{u—), b = riij{u—) = 3Cf{u—)Cj{u—) + 3Ci{u—)Cj{u—) 

and noting that a < 21^, b < 6/^ and 

Z{u) = Y^{u)s3{u) < Y\u)In(u), (7.21) 

we have 

,-- , 6a^Z{u-)Y^{u-) 3abY^{u-) 6aHY^{u-) 30aY^{u-)I''{u-) 72aY^{u-)I^{u-) 

\Ra,bV^-)\ < 9 1 1 o < o 1 o ■ 

With these observations we can represent Z in terms of stochastic integrals with respect to Ve as fol- 
lows. Recall 7^1,7^2*'"'^ introduced in the proof of Lemma [7.2[ Define ai{u) = C(2,6,ii) and a2 {u) = 
C{9ij{u),rjij{u),u). Also, let 

Ueiu) = di(u)l^^(„)(e), Ul'^iu) = d2^(u)l^(,,,)^^^(e). 
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Then 




i<j 



(7.22) 



Recalling that Ve is a rate Poisson process, one can write Z as 




eG£2 \ i<j 

Also, once again using independence of Poisson processes Ve, 




L 



B{s)ds. 



The proof is now completed upon using ()7.10p and (17. lip as for Lemma 17.21 ■ 

As is clear from the above lemma, a precise analysis of Z will involve considering several terms of the 
form I^Y^. The following lemma shows that such terms are asymptotically negligible for suitable 9,^. 

Lemma 7.8 For any 0, > and p > satisfying j{26 — i) — 1) < p 



as n ^ CO. 

Proof: From Proposition 15. H 

{In{t) < B{logn)^/{tc - tf for ah t < t„} occurs whp. 
Also, by Proposition 17.11 and (j7.4p we know that 

sup|y(t)-y(t)| =Op(y(t„)) 

t<t„ 

and from ()7.2p . for t near tc, y{tc) ~ {tc — t)/a. Using these bounds in the above integral proves the 
result. ■ 

Proof of Proposition [776l By the integral representation of the process Z given in Lemma [7. 71 and 
the fact that z solves the differential equation ()7.17p . we have that 





Jo Jo 




(7.23) 



Now the integrand in the first term can be bounded as 



3\x^{u)Y^{u) - x2(n)y^(n)| < SY'"^ {u)\x^ (u) - x'^{u)\ + 3x^{u)\Y^{u) - /(u)| 

< 6|x(n) - x{u)\ + 9\Y{u) - y{u)\. 



(7.24) 
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The integrand in the second term in (j7.23p can be decomposed as 

z{u)y{u)x^ {u) — Z{u)Y{u)x'^{u) = {u)y{u){z{u) — Z{u)) 

+ y{u)Z{u){x^{u) - x2(n)) + Z{u)x'^{u){y{u) - Y{u)). 

Thus the second integral in ()7.23p can be bounded by 



\{z{u) - Z{u))\du + Q j Z{u)\x{u) - x{u)\du + 2> I Z{u)\Y{u) - y{u)\du. (7.25) 
Combining (j7?25|) and ([7^ we get that 

\Z{t) - z{t)\ < e„ + 3 / \{Z{u) - z{u))\du 

where 

En = 9tc sup (|x(s) - x{s)\ + \Y {s) - y{s)\) + 6 / Z{u)\x{u) - x{u)\du + 3 / Z{u)\Y{u) - y{u)\du 

s<t„ 



pin pin 

+ / \Ri{u)\du+ \R2{u)\du + sup \M{t)\ 

Jo Jo i<in 



= m + m + vs + v-i + + m- 

By Gronwall's lemma, it is enough to show that e„ — )• in probability as n — )• oo. Let us show each of the 
six constituents of e„ satisfy this asymptotics. By Lemma 16.41 and Proposition 17. II 771 —t- in probability. 
Again by Lemma 16.41 for any {} < 1/2, whp, 

m<—si Z{u)du < — o / Y'^{u)In{u)du. 
n Jo ^ Jo 

Using Lemma EU the last term converges to in probability as n — )• cxd. Thus 7/2 —t- in probability. 

An identical argument, using Proposition 17. II instead of Lemma 16.41 shows that T/3 — )■ in probability. 
For r/4, note that from (17.19p . 

n 

Lemma 17.81 now shows that r/4 — )• in probability. A similar argument, using the bounds in Lemma 17.71 
on R2{u) establishes that 7/5 — )• in probability. For rjQ note that for an arbitrary stopping time r 

tAr 



pt/\T 

E sup |Mt|^ < 4E[(M)(t At)] =4E / B{u)du. 

i<inf\T Jo 



The bound on B{u) in (I7.20p along with Lemma 17.81 and a localization argument similar to the one used 
in the proof of Lemma 17.41 now shows that ryg converges to in probability. 

7.4 Proof of Proposition 14.11 

We now complete the proof of Proposition 14.11 Proof of (j4.2p follows from Proposition 17.11 and the 
discussion immediately above the proposition. Proof of ()4.ip is immediate from Proposition 17.61 Finally 
we consider (14. 3p . From Proposition 17.11 and (17. 4p 



an 



1+7 
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Also, from Proposition 15.11 

Vn^T log'^n 

as n — 7- oo. Combining, and recalling that 7 G (1/6, 1/5), we have 

52 (tn) 

This completes the proof of Propostion 14.11 ■ 

8 Proof of Theorem EOl 

We will now complete the proof of Theorem 12.11 As always, we write the component sizes as 

and write the scaled component sizes as 

^n"W = {^,Cn [tc + ^'"o^^^ ■■^>l^^ (C.(A) : i > 1) (8.1) 
Then Proposition 14.11 proves that with 

^"7+1/3 

and 7 S (1/6, 1/5) we have, as n — )• 00, 

Z^i (Cj(An)) p 1 I \ f n Ci(A„) 



EaC.(A„))^]' E.(C.(A„))' EaC.(A„))^ 

We shall now give an idea of the proof of the main result, and postpone precise arguments to the next 
two sections. The first step is to observe that the asymptotics in (18. 2p imply that the C^f process at time 
An satisfies the regularity conditions of Proposition 4 of t4j . The second key observation is that the scaled 
components merge in the critical window at a rate close to that for the multiplicative coalescent. Indeed, 
note that for any given time t components i < j £ BF(t) merge in a small time interval [t, t + dt) at rate 

-{l-x\t))Ci{t)C,{t). 
n 

Thus letting A = (t — tc)n^^^ / (a/3^/^) be the scaled time parameter, in the time interval [A, A + dX), these 
two components merge at rate 



a 



l-x2(t, + /32/3a-A_) )c,(A)C,(A). 



Now since, for large n, 



n 
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and from [26], a(l - x^{tc)) = 1 (see ([23])) we get 



7y(A) «C,(A)C,(A) 

which is exactly the rate of merger for the multiphcative coalescent. The above two facts allow us to 
complete the proof using ideas similar to those in [7]. Let us now make these statements precise. 

As before, throughout this section tn = tc — n~'^ = tc + fi'^^^oi^^j^, where 7 is fixed in (1/6, 1/5). We 

will first show that C'^^(A) — ^ ^{^) i^i ^\ each A G M and at the end of the section show that, in 

fact, C"^^ — ^ X in D((— cxd,oo) : /^). Now fix A G M. By choosing n large enough we can ensure that 
A > A„. Henceforth consider only such n. Recall that COMn{t) denotes the subgraph of BF„(t) obtained 
by deleting all the singletons. Let Yl,i(^coM denote the summation over all components in COAin, and 
denote the summation over all components in BF„. Since 

Y.mX)f- E m>^)f<^,Y.^ = 0{Vn"% (8.3) 

i ieCOM i=l 

it suffices to prove Theorem 12.11 and verify Proposition 14.11 with BF„(t) replaced by COAinit)- We write 
instead of "^i^coM ^'-'^ simplicity of the notation from now on. We begin in Section ISTl by constructing 
a coupling of {COM.n(t)}t>tn with two other random graph processes, sandwiching our process between 
these two processed, and proving statements analogous to those in Theorem 12.11 for scaled component 
vectors associated with these processes. Proof of Theorem 12.11 will then be completed in Section [8^ 



8.1 Coupling with the multiphcative coalescent 

Lower bound coupling: Let, for t > t„, CC'7W^(t) be a modification of CC'A^n(i) such that CC'A4^(t„) = 
COAinitn), and when t > t„, we change the dynamics of the random graph to the Erdos-Renyi type. 
More precisely, recall from Section [5.21 that a jump in BF„(t) can be produced by three different kinds of 
events. These are described in items (i), (ii) and (iii) in Section [5]2j COAi~{t), t > tn is constructed from 
COAi~{tn) by erasing events of type (i) and (ii) (i.e. immigrating doubletons and attaching singletons) 
and changing the probability of edge formation between two non-singletons (from that given in (|5.10p ) to 
the fixed value 5* (t„)/n. Since 6* (t) is nondecreasing in t, we have that COMn{tn + •) COA4~{tn + •)• 
Denote by (A) = (Cj~(A) : i > I) the scaled (as in (18. Ij) ) component size vector for CO^A~{t). From 
Proposition 4 of [4J, it follows that for any A G M, 

C;(A)AX(A) (8.4) 

in Indeed, note that the first and third convergence statements in ()8.2p hold with Cj replaced with 
C~ since the contributions made by singletons to the scaled sum of squares is 0{n~^/^) (see ()8.3p ) and to 
the sum of cubes is even smaller. This shows that the first and third requirements in Proposition 4 of [1] 
(see equations (8), (10) therein) are met. To show the second requirement in Proposition 4 of [3], using 
the second convergence in ()8.2p . 



lim (n^/^p-y^Y mAP^lMX -Xn) I \ 

n^oo\\ ) n nV3 E.(CT(An))V 

lim ab*^{tn)X - XnioKitn) - 1) 

n— >oo 

A - lim An(a6*(tn) - 1), 
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where the last equahty follows on observing that, as n — )• oo, Kitn) ^ 1 - a;2(tj and a{l - ^^(tj) = 1. 
Also, 

^-7+1/3 

lim An|a6* (fn) - 1| = lim — -— — 16* - a | 

^-7+1/3 



^2/3 l^o(^fa)) ~ ^o(a;(tc))| 



< di lim n-T+^/3|x(t„) -a;(tc)| 

< hm d2 - + n-^+^/^\tn - t^) , 

where the second equality follows from (j5.12p . The first term on the last line converges to using Lemma 
16.41 For the second term note that n~^^'^l'^\tn — ic| = n~'^^'^l'^n~'^ which converges to since 7 > 1/6. 
Thus we have shown that the expression in ()8.5p converges to A as n — )• 00 and therefore the second re- 
quirement in Proposition 4 of [1] (see equation 9 therein) is met as well. This proves that C~{\) X{\) 
in for every A G M. Although Proposition 4 of [3] only proves convergence at any fixed point A, from 
the Feller property of the multiplicative coalescent process proved in Proposition 6 of the same paper it 
now follows that, in fact, — ^ X. in 00,00) : /|). 



Upper bound coupling: Let us construct {COA^^(t) : t > t„} in the following way. Let = 

tc + n^'^ and let 

A+ = (4 - tc)nV3/(a/32/3) = nV3-7/(«^2/3)_ 

Let COAi^{tn) be the graph obtained by including all immigrating doubleton and attachments during 
time t G [tn^^n] to graph ol COAin{t), along with all the attachment edges. Namely, we construct 
COM.^{tn) by including in COMn{tn) all events of type (i) and (ii) of Section [512] that occur over [tn, t'^]- 
For t > tn the graph evolves in the Erdos-Renyi way such that edges are added between each pair of vertices 
in the fixed rate 6*(t+)/n. The coupling between COAi^{- + tn) and COA^„(- + t„) can be achieved as 
follows: Construct a realization of {COM.n{t) '■ tn <t < t+} first, then use hn{t'^) — hn{t) to make up for 
ah the additional edges in CC'X+(t) for tn < t < t+. Note that COMn{tn + •) <d COM^{tn + •) over 
[0,4 -in]. 

Let Cn{\) = (C+(A) : i > 1) be the scaled (as in (|8.ip ) component size vector for COAin- We will once 
more apply Proposition 4 of We first show that the three convergence statements in (18. 2p hold with 
Cj replaced with Cf. For this it will be convenient to consider processes under the original time scale. 
Write Cn\tn) = €{. Also denote by {C^} the component vector obtained by adding all events of type (ii) 
only, to COAinitn) (i-e. attachment of singletons to components in COAinitn)), 

over Since c* is 

bounded by 1, is stochastically dominated by the sum of Cj independent copies of Geometric(p), with 
p = e*"-*" = e~^" ^. Thus 

Ui cf — Ci <d Negative-binomial (r,p) with r = Ci,p = e~^" ~' . 

The random graph COA4f{tn) contains components other than {Cf}- These additional components 
correspond to the ones obtained from doubletons immigrating over [tn , ] • Since there are at most ti 
vertices, the number of such doubletons is bounded by n/2. Denote by {Cf}^^ the components 
corresponding to such doubletons. Once again using the fact that c* < 1, we have that 

C+ <d 2 Negative-binomial (2, p) with p = e"^"~^. 

Write 

AT 

Sk = ^(C,)^ 5+ = ^(C+)'= + ^(C+)'' for A; = 2, 3 and / = maxC,, /+ = max{maxC+, maxC+}. 

i i i=X 
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The following proposition shows that Propostion 14.11 holds with {S2{tn),Ss{tn),Cn\tn)) replaced with 

{S+{tn),S+{tn),I+{tn)). 



Proposition 8.1 As n — )• oo, 



n 



SI 
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/+ = 6(1) 
^ 1 



1 



4/3 ( ^ _ 1 ] ^ 0. 



Proof: An elementary calculation shows that if U is Negative-binomial (r, e ^" ^) then for some di G 
(0,00) 

p(c/ > 37-V) < % 

and thus, as n — )• 00, 

P(maxC+ > (1 + 37"^)/) < ¥{ui > 37"^Ci for some i = 1, • • • n) ^ 0. 

i 

A similar calculation shows that, for some ^2 € (0,oo), as n — )• 00. 

P( max Ct > d2) 0. 
\=l,...Ar * 

The first statement in the proposition now follows on combining the above two displays. 
Next, note that for Negative-binomial(r,p), the first, second and third moments are 

Ml = -r(l - p) 
P 
1 



M2 = ^[r'^{l-pf +r{l-p)] 



p. 
1 



M3 = ^[r3(l + 3r2(l + r(4-9p + 7y>2 - 2/)]. 

Prom ([7:26]1 and gU it follows that ^2 = Q{n^+'^) and S3 = G(n^+3'^). Also, clearly, Ylid = 0{n). 
Write D2 = S+ -S2 = E^=l(C+)' + Eii^CiU, + nf), then 

nD2\{C,h]<d2 (n-n-^ + Y,[{Cifn-^ + {Cifn-^^ +C,n-'']j = 0{n) 

thus D2/S2 and consequently S2 /S2 1- 

Write D3 = S^ -S3 = EZii^tf + Ei[3(C»)2u, + 3C,n2 + nf]. One can similarly show that 



E[D3mU = 0{n 

thus D3/S3 and so 5^/53 1. 

To prove the third convergence, it suffices to prove 



i+2r\ 



^-0. (8.6) 
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By the asymptotics shown above, we have 

= 0(n4/3+i-2{i+7)) = 0(ni/3-27) 

As 7 > 1/6, (18. 6p follows and thus the proof is completed. ■ 

For scaled component size vector of COAi^, the above proposition shows that the statements in (j8.2p 
hold with Ci replaced with . In particular, the first and third requirements in Proposition 4 of [4J are 
met by {C^} Also, using the second convergence in (j8.2p . a calculation similar to that for (j8.5p shows 
that 

n^coy\ I n ni/s Ei(CT(^»)) / 

Therefore the second requirement in Proposition 4 of ^ is satisfied. This proves that 



Cl(\)^X(\). (8.7) 



in for every A G M. Using Proposition 6 of [4J once again it now follows that — )■ X in !D((— oo, cxd) : 



8.2 Completing the proof of Theorem 12.11 

By [6l[7|, there is a natural partial order -< on Informally, interpreting an element of /| as a sequence of 
cluster sizes, x, y G x ^ y if y can be obtained from x by adding new clusters and coalescing together 
clusters. The coupling constructed in Section [8T] gives that, for every, A G (A„,A+) 

Since, as n — )• cx). An — >• -co and A^ — t- +oo, ()8.4p . ()8.7p along with Lemma 15 of [7J yield that 

ar(A) AX(A) 

for all A G M. 

Finally we argue convergence in P((— cxd,C)o) : Z^). For x,y G let d^(x, y) = ~ Vif' ■, 

X = {xj}, y = {yi\. Then d^(x,y) < Y^iuf — Yli^l whenever x ^ y. To prove that — )• X in 



cxD,oo) : if) it suffices to prove that 



sup d(C^^, C„ ) 0, for all — oo < Ai < A2 < 00. 



AG[Ai,A2] 

Fix Ai,A2 as above. Then 

sup d(Cr,C;)< sup E(C+(A))2-5](Cr(A))^]. (8.9) 

AG[Ai,A2] Ae[Ai,A2] ^ ~^ 

Let, for A G M, 

U+{X) = ^(C7(A))2, U.{X) = Y.{C.{X)f and V{X)=U+{X) -U.{X). 

i i 

From Lemma 15 of [7J, V(A) for every A G M. Thus it suffices to show that V is tight in 'D((— 00, 00) : 
M+). Note that both and U- are tight in D{{— 00,00) : ^+)- Although, in general difference of 
relatively compact sequences in the P-space need not be relatively compact, in the current setting due to 
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properties of the multiplicative coalescent this difficulty does not arise. Indeed, if {X^{t),t > 0} denotes 
the multiplicative coalescent on the positive real line with initial condition x E then, for 5 sufficiently 



small 



sup E(d2(X^(r),x) Al) <E 

rgr((S) 



< 2 > 6xiXj ■ 2xiXj < 25||x| 



i<j 



where, ||x|| = xf)^^'^ , T{5) is the family of all stopping times (with the natural filtration) bounded 
by 5. Using the above property, the Markov property of the coalescent process and the tightness of 
sup;^g[;^^ ^+(-^)) ^^V\^[\i,\2\^-{^) ^au Verify Aldous' tightness criteria (see Theorem VI. 4. 5 in [24J) 
for V thus proving the desired tightness. ■ 
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